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Abstract We present the program package GoSam 
which is designed for the automated calculation of one- 
loop amplitudes for multi-particle processes in renorma- 
lisablc quantum field theories. The amplitudes, which 
are generated in terms of Feynman diagrams, can be 
reduced using cither D-dimcnsional integrand-level de- 
composition or tensor reduction. GoSam can be used to 
calculate one-loop QCD and/or electroweak corrections 
to Standard Model processes and offers the flexibility 
to link model files for theories Beyond the Standard 
Model. A standard interface to programs calculating 
real radiation is also implemented. We demonstrate the 
flexibility of the program by presenting examples of pro- 
cesses with up to six external legs attached to the loop. 

Keywords NLO calculations ■ automation ■ hadron 
colliders 
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1 Introduction 

The Standard Model is currently being re-discovered 
at the LHC, and new exclusion limits on Beyond the 
Standard Model particles - and on the Higgs mass - 
are being delivered by the experimental collaborations 
with an impressive speed. Higher order corrections play 
an important role in obtaining bounds on the Higgs bo- 
son and New Physics. In particular, the exclusion limits 
for the Higgs boson would look very different if we only 
had leading order tools at hand. Further, it will be very 
important to have precise theory predictions to con- 
strain model parameters once a signal of New Physics 
has been established. Therefore it is of major impor- 
tance to provide tools for next-to-leading order (NLO) 
predictions which are largely automated, such that sig- 
nal and background rates for a multitude of processes 
can be estimated reliably. 

The need for an automation of NLO calculations 
has been noticed some time ago and lead to public pro- 
grams like FeynArts Q] and Q Graf [5] for diagram gen- 
eration and FormCalc/LoopTools[3J and GRACE gj 
for the automated calculation of NLO corrections, pri- 
marily in the electroweak sector. However, the calcu- 
lation of one-loop amplitudes with more than four ex- 
ternal legs were still tedious case-by-case calculations. 
Only very recently, conceptual and technical advances 
in multi-leg one-loop calculations allowed the calcu- 
lation of six-point [^ElLTlIHlElIinillllllSlIISlIIllIISlIinWIZl 
[T51[TW^I^2P5P11 and even seven-point gMM Pro- 
cesses at all, and opened the door to the possibility 
of an automated generation and evaluation of multi-leg 
one-loop amplitudes. As a consequence, already exist- 
ing excellent public tools, each containing a collection 
of hard-coded individual processes, like e.g. MCFM [27J 
US], VBFNLO EMU], MC@NLO [32132], POWHEG-Box 
[3JI3J], POWHEL [311137], can be flanked by flexible 
automated tools such that basically any process which 
may turn out to be important for the comparison of 
LHC findings to theory can be evaluated at NLO accu- 
racy. 

We have recently experienced major advances in 
the activity of constructing packages for fully auto- 
mated one-loop calculations, see e.g. [38, 39 ,40 ,4111421 
|4"3] . The concepts that lead to these advances have 
been recently reviewed in [44] . Among the most im- 
portant developments are the integrand-reduction tech- 
nique [451146] and the generalized n-dimensional unitar- 
ity [47] . Their main outcome is a numerical reconstruc- 
tion of a representation of the tensor structure of any 
one-loop integrand where the multi-particle pole con- 
figuration is manifest. As a consequence, decomposing 
one-loop amplitudes in terms of basic integrals becomes 



equivalent to reconstructing the polynomial forms of 
the residues to all multi-particle cuts. Within this al- 
gorithm, the integrand of a given scattering amplitude, 
carrying complete and explicit information on the cho- 
sen dimcnsional-rcgularisation scheme, is the only in- 
put required to accomplish the task of its evaluation. 
In fact, the integration is substituted by a much simpler 
operation, namely by polynomial fitting, which requires 
the sampling of the integrand on the solutions of gen- 
eralised on-shell conditions. 

In this article, we present the program package Go- 
Sam which allows the automated calculation of one- 
loop amplitudes for multi-particle processes. Amplitudes 
are expressed in terms of Feynman diagrams, where 
the integrand is generated analytically using QGRAF [5J, 
FORM g3J, spinney g5] and haggies [50]. The individ- 
ual program tasks are steered via python scripts, while 
the user only needs to edit an "input card" to specify 
the details of the process to be calculated, and launch 
the generation of the source code and its compilation, 
without having to worry about internal details of the 
code generation. 

The program offers the option to use different reduc- 
tion techniques: either the unitarity-based integrand re- 
duction as implemented in Samurai gD] or traditional 
tensor reduction as implemented in Golem95C [5T1I52] 
interfaced through tensorial reconstruction at the in- 
tegrand level [S3] , or a combination of both. It can be 
used to calculate one-loop corrections within both QCD 
and electroweak theory. Beyond the Standard Model 
theories can be interfaced using FeynRules [54] or Lan- 
HEP [55] . The Binoth-Les Houches- interface [56] to pro- 
grams providing the real radiation contributions is also 
included. 

The advantage of generating analytic expressions for 
the integrand of each diagram gives the user the flexi- 
bility to organize the computation according to his own 
efficiency preferences. For instance, the computing al- 
gorithm can proceed either diagram-by-diagram or by 
grouping diagrams that share a common set of denom- 
inators (suitable for a unitarity-based reduction), and 
it can deal with the evaluation of the rational terms ei- 
ther on the same footing as the rest of the amplitude, or 
through an independent routine which evaluates them 
analytically. These options and the other features of 
GoSam will be discussed in detail in the following. 

In Section [5] after giving an overview on the di- 
agram generation and on processing gauge-group and 
Lorentz algebra, we discuss the code generation and 
the reduction strategies. The installation requirements 
are given in Section [3J while Section 2] describes the us- 
age of GoSam, containing all the set-up options which 
can be activated by editing the input card. In Section 
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[5] we show results for processes of various complexity. 
The release of GoSam is accompanied by the generated 
code for some example processes, listed in Appendix Kl 



2 Overview and Algorithms 

2.1 Overview 

GoSam produces, in a fully automated way, all the code 
required to perform the calculation of one-loop matrix 
elements. There are three main steps in the process of 
constructing the code: the generation of all contributing 
diagrams within a process directory, the generation of 
the Fortran code, and finally compiling and linking 
the generated code. These steps are self-contained in 
the sense that after each step all the files contained in 
the process directory could be transfered to a different 
machine where the next step will be carried out. 

In the following sections we focus on the algorithms 
that are employed for the construction of the code to 
produce and evaluate matrix elements. 

The first step (setting up a process directory), which 
consists in the generation of some general source files 
and the generation of the diagrams, is described in Sec- 
tion [521 The second step (generating the fortran code) 
is carried out by means of advanced algorithms for al- 
gebraic manipulation and code optimization which are 
presented in Sections 12.31 and 12.41 The third step (com- 
pilation and linking) is not specific to our code genera- 
tion, therefore will not be described here. 

The practical procedures to be followed by the user 
in generating the code will be given in Section 21 which 
can be considered a short version of the user manual. 



2.2 Generation and Organisation of the Diagrams 

For the diagram generation both at tree level and one- 
loop level we employ the program QGRAF [2] . This pro- 
gram already offers several ways of excluding unwanted 
diagrams, for example by requesting a certain number 
of propagators or vertices of a certain type or by speci- 
fying topological properties such as the presence of tad- 
poles or on-shell propagators. Although QGRAF is a very 
reliable and fast generator, we extend its possibilities 
by adding another level of analysing and filtering over 
diagrams by means of Python. This gives several advan- 
tages: first of all, the possibilities offered by QGRAF are 
not always sufficient to distinguish certain classes of di- 
agrams (see examples in Fig.[IJ; secondly, QGRAF cannot 
handle the sign for diagrams with Majorana fcrmions 
in a reliable way; finally, in order to fully optimize the 




u(k 2 ) 



u(k>2 

(a) Diagram 1. 




(b) Diagram 2. 



Fig. 1 Two examples for diagrams which are difficult to iso- 
late using QGRAF. The diagram in Fig. |l(a)| is zero in dimen- 
sional regularisation. However, in QGRAF there is no operator 
to identify this type of diagrams. In Fig. 1 1(b)] the Z boson is 
emitted from a closed quark line. These diagrams form a sep- 
arate gauge invariant class and could be treated separately 
from diagrams where the Z boson comes from an external 
quark line. 



reduction, we want to classify and group diagrams ac- 
cording to the sets of their propagators. 

Within our framework, QGRAF generates three sets 
of output files: an expression for each diagram to be 
processed with FORM [48] . Python code for drawing all 
diagrams, and Python code for computing the proper- 
ties of each diagram. The information about the model 
for QGRAF is either read from the built-in Standard Mo- 
del file or is generated from a user defined LanHEP [55] 
or Universal FeynRules Output (UFO) [51] file. 

The Python program automatically performs sev- 
eral operations: 

— diagrams whose color factor turns out to be zero are 
dropped automatically; 

— the fermion flow is determined and used to compute 
an overall sign for each diagram, which is relevant 
in the presence of Majorana fcrmions; 

— the number of propagators containing the loop mo- 
mentum, i.e. the loop size of the diagram, the tensor 
rank and the kinematic invariants of the associated 
loop integral are computed; 

— diagrams with an associated vanishing loop integral 



(see Fig. 1(a) I are detected and flagged for the dia- 
gram selection; 
— all propagators and vertices are classified for the di- 
agram selection; diagrams containing massive quark 
self-energy insertions or closed massless quark loops 
are specially flagged. 

Any one-loop diagram can be written in the form 
d n q Af(q) 



V 



T n/2 



Uli [(q + rty-mf + iS] 
where the numerator is a polynomial of tensoiQ rank 



(1) 



Af(q) = Co + C? %1 +... + C^-^q^ ■ ■ ■ q, r 



(2) 



1 Index contractions in Eq. ([2J are understood in n- 
dimensional space. 
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and the N x N kinematic matrix is defined as 

Sij = (n - r 3 f - mf - uij . (3) 

All masses can be either real or complex. Important 
information about the integrals that will appear in the 
reduction of each one-loop diagram is contained in the 
tensor rank r of the loop integral and its kinematic 
matrix Sij. 

We define a preorder relation on one-loop diagrams, 
such that T>i ^ T> 2 if their associated matrices S(T>i) 
and S(T>2) are related by a finite (not necessarily unique) 
chain of transformations 

S(X> 2 ) AS'A..A S<Pi) (4) 
where each transformation is one of the following: 

— the identity, 

— the simultaneous permutation of rows and columns, 

— the simultaneous deletion of the row and column 
with the same index, which corresponds to pinching 
the corresponding propagator in the diagram. 

The relation "X" can be read as "appears in the reduc- 
tion of" . Our algorithm groups the one-loop diagrams 
T>i, . . . , T>d of a process into subsets V\, . . . , Vq such 
that 

— V\, . . . , Vq form a partition of {2?i, . . • and 

— each cell Vi contains a maximum element max Vi £ 
Vi, such that V ^ maxVi ,VX> e V { . 

The partitioning procedure provides an important gain 
in efficiency, because while carrying out the tensor re- 
duction for the diagram maxl^, all other diagrams in 
the same cell Vi are reduced with virtually no addi- 
tional computational cost. The gain in efficiency can 
be observed when reducing the diagram using the OPP 
method [55J and its implementations in Cut Tools [57] 
and Samurai [JO], as well as in classical tensor reduc- 
tion methods as implemented e.g. in Golem95C [5T1I52] . 
PJFRY [5S] and LoopTools [3J[5J]. 

In order to draw the diagrams, we first compute an 
ordering of the external legs which allows for a planar 
embedding of the graph. Such ordering can always be 
found for a tree or a one-loop graph since non-planar 
graphs only start to appear in diagrams with two or 
more loops. After the legs have been assigned to the 
vertices of a regular polygon, we use our own implemen- 
tation of the algorithms described in [SD] for fixing the 
coordinates of the remaining vertices; the algorithm has 
been extended to determine an appealing layout also 
for graphs containing tadpoles. Starting from these co- 
ordinates and using the package Axodraw [HT], GoSam 
generates a I^TjrjX file that contains graphical represen- 
tations of all diagrams. 



2.3 Algebraic Processing 
2.3.1 Color Algebra 

In the models used by GoSam, we allow one unbroken 
gauge group SU(Nc) to be treated implicitly; any ad- 
ditional gauge group, broken or unbroken, needs to be 
expanded explicitly. Any particle of the model may be 
charged under the SU (Nc) group in the trivial, (anti- 
fundamental or adjoint representation. Other represen- 
tations arc currently not implemented. 

For a given process we project each Feynman dia- 
gram onto a color basis consisting of strings of gener- 
ators T Al T A f ■■■T Ap . and Kronecker deltas Sa but 
no contractions of adjoint indices and no structure con- 
stants f ABC . Considering, for example, the process 

u(l) + u(2)->Z(3)+ff(4) + <7(5) 

GoSam finds the color basis 

= q^q^g^g^MT^T^}, 

where q\'^ and g A t ' are the color parts of the quark 
and gluon wave functions respectively. The dimension 
of this color basis for N g external gluons and N q q quark- 
antiquark pairs is given by [55] : 

d(N g , N qq ) = ( 7 ) • (N g + N qq - i)\ . (5) 

i=o ^ ' 

It should be noted that the color basis constructed in 
this way is not a basis in the mathematical sense, as one 
can find linear relations between the vectors |c,} once 
the number of external partons is large enough. 

Any Feynman diagram can be reduced to the form 

k 

V = Y / C l \c l ) (6) 

i=l 

for the process specific color basis |ci), . . . , \ck) by ap- 
plying the following set of relations: 

1ij T H = T R (saSkj - -^-<%4;^ , (7) 

fABC _ jrpArpBrpC rpArpC rpB\ I q\ 

J - { 1 ij 1 jk 1 ki '.;';/•'/.'• \ S ) 

The same set of simplifications is used to compute the 
matrices (ci\cj) and (ci\Ti - Tj\cj) . The former is needed 
for squaring the matrix element, whereas the latter is 
used to provide color correlated Born matrix elements 
which we use for checking the IR poles of the virtual 
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amplitude and also to provide the relevant information 
for paxton showers like POWHEG [fJ3"H33U34j . For the above 
example, GoSam obtains^] 



(ci\cj) = T R C F 



\N C -1) -1 



N c 

-1 (JVg - 1) N c 
N c N c N C/ 



(9) 



Similarly, the program computes the matrices (cj|Tj • 
Tj\cj) for all pairs of partons I and J. 

If jM ^ denotes the tree- level matrix element of the 
process and we have 



(10) 



then the square of the tree level amplitude can be writ- 
ten as 



mm = £ (c| 0) Vd 0) 



Ci|Cj). 



(11) 



For the interference term between leading and ncxt- 
to-leading order we use a slightly different philosophy. 
First of all we note that it is sufficient to focus on a 
single group V a as defined in Section l2~2l 



E /— 



d™g AA a (g) 



i/2 ttW 



nii[(9 + n) 



iS] 



+ h.c. (12) 



In order to reduce the complexity at the level of the 
reduction, we perform the contraction with the tree- 
level already at the integrand level, 



In GoSam we have implemented the 't Hooft-Velt- 
man scheme (HV) and dimensional reduction (DRED). 
In both schemes all external vectors (momenta and po- 
larisation vectors) are kept in four dimensions. Internal 
vectors, however, are kept in the n-dimcnsional vector 
space. We adopt the conventions used in [15], where k 
denotes the four-dimensional projection of an in gen- 
eral n-dimensional vector k. The (n — 4)-dimensional 
orthogonal projection is denoted as k. For the integra- 
tion momentum q we introduce in addition the symbol 
fi 2 = — q 2 , such that 



2 -2,-2 -1 2 

q =q +q =q - Ijl . 



(14) 



We also introduce suitable projectors by splitting the 
metric tensor 



n fi = 0, ft = 4, ft = n-4. 

(15) 



In the follwing, we describe the 't Hooft algebra 
in detail. For DRED, the only differences are that the 
numerator algebra is performed in four dimensions for 
both external and internal vectors (i.e. q = q) and that 
in the very end all appearances of q 2 are replaced by 

-2 2 

T — M ■ 

Wave Functions and Propagators GoSam contains a 
library of representations of wave functions and propa- 
gators up to spin twcH. The exact form of the interac- 
tion vertices is taken from the model files. 

The representation of all wave functions with non- 
trivial spin is based on massless spinors. Each massive 
external vector pi is replaced by its light-cone projec- 
tion li with respect to a lightlike reference vector k, 



Af a (q)=-£ ( Cl | Cj )(cfyC<%), 

where (jj is formed by the sum over the corresponding 
coefficients of all diagrams T> E V a . 



(13) 



2.3.2 Lorentz Algebra 

In this Section we discuss the algorithms used by Go- 
Sam to transform the coefficients and C t - ((?), as 
defined in the previous section, such that the result 
is suitable for efficient numerical evaluation. One of 
the major goals is to split the n-dimcnsional algebra 
(n = 4 — 2e) into strictly four-dimensional objects and 
symbols representing the higher-dimensional remain- 
der. 



2 In the actual code the results are given in terms of Tr 
and Nq only. 



2pi ■ k 



(16) 



For spin 1 /2 particles we use the assignment of wave 
functions as shown in Table [T] here, we quote the def- 
inition of the massive spinors from |49] assuming the 
splitting of Eq. (fl6]l: 



|p ± > = IO± 

(p ± \ = (l\± 



[Ik] 



[kl] 



± 



± 



(ik) 



(ki) 



\k), (17a) 



(k\. (17b) 



In order to preserve the condition that for any loop inte- 
gral the tensor rank does not exceed the number of loop 



3 Processes with particles of spin 3 / 2 an d spin 2 have not been 
tested extensively. Furthermore, these processes can lead to 
integrals where the rank is higher than the loop size, which 
at the moment are neither implemented in Samurai nor in 
Golem95C. 
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(a) Assignment of initial and fi- 
nal states for quarks and leptons. 





I , q 


l+,q 


initial 




Va(k,j 3 ) 


final 







(b) Wave functions for mass- 
less fermions. 
u a (k,+l) = v a (k,~l) = \k) 
u a (k,-l) = v a {k,+l) = \k] 
u a {k,+l) = v a {k,-l) = [k\ 
u a (k,-l) = v a {k,+l) = (k\ 



(c) Wave functions for massive fermions. 



u a (p, +1) = 
u a (p, -1) = 

V a (p, +1) = 
V a (p, -1) = 



U a (p, +1) = [p^ 

■a a (p, -l) = (p~ 

Va{p, +1) = (p~ 
Va(p, -1) = [p~ 



Table 1 Assignment of quark and lepton wave functions. We label the physical spin states by jz = ±1, which is twice the 
3-component of the spin. The wave functions assigned in Table (a) are mapped onto the bracket notation used in spinney | 49 | 
as defined in Tables (b) and (c). 



propagators we fix all gauge boson propagators to be in 
Feynman gauge. Their wave functions are constructed 
as El 



(g|7> b ] 



V2 [p b q] ' 



(18) 



where p b = p in the massless case and p = I according 
to Eq. (JTHJ) in the massive case. In the latter case the 
third polarisation is defined as 

e^p i 0) = J-(2p\-p,\. (19) 



The wave functions and propagators for spin 3 /2 and 
spin 2 particles correspond to those in [55] . 

Simplifications Once all wave functions and propaga- 
tors have been substituted by the above definitions and 
all vertices have been replaced by their correspond- 
ing expressions from the model file then all vector-like 
quantities and all metric tensors are split into their 
four-dimensional and their orthogonal part. As we use 
the 't Hooft algebra, 75 is defined as a purely four- 
dimensional object, 75 = ie AII yp (T 7 M 7 iy 7 p 7' 7 . By applying 
the usual anti-commutation relations for Dirac matri- 
ces we can always separate the four-dimensional and 
(n — 4)-dimcnsional parts of Dirac traces, as we can use 
the fact that I62ll49l 



tr(l) -tr^j ■■■7 w 7/i i+ i ■■'Tm+J = 

tr(7 Pl -tr(7 Mj+1 •••%,+„)• (20) 

The same logic applies to open spinor lines such as [49] 

tr (!) ' •••%,%,+! •••7w+pl fc 2> = 

■••7wl fc 2) -tr(> !+1 ■■■% l+p ). (21) 

While the (n — 4)-dimensional traces are reduced com- 
pletely to products of (n — 4)-dimensional metric ten- 
sors g^ v , the four-dimensional part is treated such that 
the number of terms in the resulting expression is kept 
as small as possible. Any spinor line or trace is bro- 
ken up at any position where a light-like vector ap- 
pears. Furthermore, Chisholm idenities are used to re- 
solve Lorcntz contractions between both Dirac traces 



and open spinor lines. If any traces remain we use the 
built-in trace algorithm of FORM gS] • 

In the final result we can always avoid the explicit 
appearance of Levi-Civita tensors, noticing that any 
such tensor is contracted with at least one light-like 
vectoiQ fc M , and we can replace 



■7 (WWMfc) - (k\%%%\k]) . (22) 



Hence, the kinematic part of the numerator, at the end 
of our simplification algorithm, is expressed entirely in 
terms of: 

— spinor products of the form (kikj), [kikj] or [kil^^lkj}- 

— dot products ki ■ kj or ki ■ q, 

— constants of the Lagrangian such as masses, widths 
and coupling constants, 

— the symbols /i 2 = q 2 — q 2 and e = (n — 4)/2. 

Treatment of R2 rational terms In our representation 
for the numerator of one-loop diagrams, terms contain- 
ing the symbols /x 2 or e can lead to a so-called R2 
term [66| , which contributes to the rational part of the 
amplitude. In general, there are two ways of splitting 
the numerator function: 



Af{q,fi 2 ,e)=No{q,H 2 )+£Nl(q,H 2 ) 



2 A/" 2 (g,/x 2 



or, alternatively, 



A/(g,/iV)= N(q)+ A7(g, M 2 ,e). 



(23a) 



(23b) 



It should be noted that in Eq. (|23a[) the terms J\fi and 
A/2 do not arise in DRED, where only terms containing 
fi 2 contribute to R2 ■ Instead of relying on the construc- 
tion of i?2 from specialized Feynman rules [rTfllrlSllMl 
[70] , we generate the R2 part along with all other contri- 
butions without the need to separate the different parts. 
For efficiency reasons, however, we provide an implicit 
and an explicit construction of the R2 terms. 



4 Any external massive vector at this point has been replaced 
by a pair of light-like ones. Contractions between two Levi- 
Civita symbols can be resolved to products of metric tensors. 
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The implicit construction uses the splitting of Eq. (|23aj) 
and treats all three numerator functions Afi on equal 
grounds. Each of the three terms is reduced separately 
in a numerical reduction and the Laurent series of the 
three results are added up taking into account the pow- 
ers of e. 

The explicit construction of i?2 is based on the as- 
sumption that each term in Af in Eq. (|23b[) contains at 
least one power of fi 2 or e. The expressions for those in- 
tegrals are relatively simple and known explicitly. Hence, 
the part of the amplitude which originates from Af is 
computed analytically whereas the purely four-dimen- 
sional part Af is passed to the numerical reduction. 



2.4 Code Generation 
2.4-1 Abbreviation System 

To prepare the numerator functions of the one-loop di- 
agrams for their numerical evaluation, we separate the 
symbol /x 2 and dot products involving the momentum 
q from all other factors. All subexpressions which do 
not depend on either q or pi 2 arc substituted by ab- 
breviation symbols, which are evaluated only once per 
phase space point. Each of the two parts is then pro- 
cessed using haggies [50], which generates optimized 
Fortran code for their numerical evaluation. For each 
diagram we generate an interface to Samurai [ID], Go- 
lem95C [5J and/or PJFRY [55]. The two latter codes 
are interfaced using tensorial reconstruction at the in- 
tegrand level [53] . 

2-4-2 Reduction Strategies 

In the implementation of GoSam, great emphasis has 
been put on maintaining flexibility with respect to the 
reduction algorithm that the user decides to use. On the 
one hand, this is important because the best choice of 
the reduction method in terms of speed and numerical 
stability can strongly depend on the specific process. 
On the other hand, we tried to keep the code flexible to 
allow further extensions to new reduction libraries, such 
that GoSam can be used as a laboratory for interfacing 
future methods with a realistic environment. 

Our standard choice for the reduction is Samurai, 
which provides a very fast and stable reduction in a 
large part of the phase space. Furthermore, Samurai 
reports to the client code if the quality of the recon- 
struction of the numerator suffices the numerical re- 
quirements (for details we refer to [10]). In GoSam we 
use this information to trigger an alternative reduction 
with either Golem95C [521 or PJFRY [55] whenever these 




Fig. 2 Reduction strategies currently implemented in Go- 
Sam: the reduction algorithm is chosen by setting the variable 
reduction_interoperation in the generated Fortran code and 
can be modified at run time. 0: Samurai only; 1: Golem95C 
only; 2: Samurai with rescue option (Golem95C); 3: Samurai 
with numerator from tensorial reconstruction; 4: same as 3 
but with rescue option(Golem95C). 11, 12 and 14 are the same 
as 1, 2, 3 (respectively) with the difference that PJFRY is used 
instead of Golem95C. 



reconstruction tests fail, as shown in Fig. [5] The reduc- 
tion algorithms implemented in these libraries extend 
to phase space regions of small Gram determinants and 
therefore cover most cases in which on-shcll methods 
cannot operate sufficiently well. This combination of 
on-shell techniques and traditional tensor reduction is 
achieved using tensorial reconstruction at the integrand 
level |53| , which also provides the possibility of running 
on-shell methods with a reconstructed numerator. In 
addition to solving the problem of numerical instabili- 
ties, in some cases this option can reduce the compu- 
tational cost of the reduction. Since the reconstructed 
numerator is typically of a form where kinematics and 
loop momentum dependence are already separated , the 
use of a reconstructed numerator tends to be faster 
than the original procedure, in particular in cases with 
a large number of legs and low rank. 

The flowchart in Fig. [2] summarizes all possible re- 
duction strategies which are currently implemented. The 
strategy in use is selected by assigning the variable 
reduction_interoperation in the generated Fortran 
code. The availability of the branches is determined 
during code generation by activating (at least one of) 
the extensions (samurai, golem95, pjf ry) in the input 
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card. Switching between active branches is possible at 
run time. In detail, the possible choices for the variable 
reduction_interoperation are the following: 

the numerators of the one-loop diagrams are reduced 

by Samurai, no rescue system is used in case the 
reconstruction test fails; 

1 the tensor coefficients of the numerators are recon- 

structed using the tensorial reconstruction at the 
integrand level, the numerator is expressed in terms 
of tensor integral form factors which are evaluated 
using Golem95C; 

2 the numerators are reduced by Samurai; whenever 

the reconstruction test fails, numerators are reduced 
using the option 1 as a backup method; 

3 tensorial reconstruction is used to compute the tensor 

coefficients; Samurai is employed for the reduction 
of the reconstructed numerator, no rescue system is 
used; 

4 as in option 3, Samurai is used to reduce the re- 

constructed numerator, Golem95C is used as backup 
option; 

11 same as 1 but PJFRY is used instead of Golem95C; 

12 same as 2 but PJFRY is used instead of Golem95C; 
14 same as 4 but PJFRY is used instead of Golem95C. 

It is difficult to make a statement about the "opti- 
mal" reduction method because this depends on the 
process under consideration. For multi-leg processes, 
e.g. bbbb production, we found that Samurai is clearly 
superior to tensor reduction in what concerns timings 
and size of the code. Concerning points which need a 
special treatment, we did not make extensive studies 
using traditional tensor reduction only, but one can 
certainly say that the combination of Samurai and 
tensorial reconstruction seems to be optimal in what 
concerns the avoidance of numerical instabilities due to 
inverse Gram determinants. 

2.5 Conventions of the Amplitudes 

In this section we briefly discuss the conventions chosen 
for the results returned by GoSam. Depending on the 
actual setup for a given process, in particular if an im- 
ported model file is used, conventions may be slightly 
different. Here we restrict the discussion to the case 
where the user wants to compute QCD corrections to 
a process and in the setup files he has put g s = 1. In 
this case, the tree-level matrix element squared can be 
written as 

\M\ 2 tIee = 4Ao = (g s ) 2b ■ a . (24) 

The fully rcnormaliscd matrix element at one-loop level, 
i.e. the interference term between tree-level and one- 



loop, can be written as 

1-^lLoop = A\A + AlA! = 2 • $t(AlA!) = 

|A4|bare + |X|ct, 5mQ + |^|ct, „. + M£f, g + |X|wf, Q = 

a s (n) (4tt) £ 2b r c_i c_ 2 , ] 

2n r(l-e) ■ {9s) ■ l Ca + T + — +0{£) \ ■ 

(25) 

A call to the subroutine Samplitude returns an array 
consisting of the four numbers (ao, Co, c_i, c_ 2 ) m this 
order. The average over initial state colours and hclici- 
ties is included in the default setup. In cases where the 
process is loop induced, i.e. the tree level amplitude is 
absent, the program returns the values for .Aj„4i where 
a factor 

( asM (4tt)* \ 2 
^ 2n r(l-e)J 

has been pulled out. 

After all UV-renormalisation contributions have been 
taken into account correctly, only IR-singularities re- 
main, which can be computed using the routine ir_sub- 
tractions. This routine returns a vector of length two, 
containing the coefficients of the single and the double 
pole, which should be equal to (c_i,c_2) and therefore 
can be used as a check of the result. 

Ultraviolet Renormalisation in QCD For UV-renormal- 
isation we use the MS scheme for the gluon and all 
massless quarks, whereas a subtraction at zero momen- 
tum is chosen for massive quarks [71] . Currently, coun- 
terterms are only provided for QCD corrections. In the 
case of electroweak corrections only unrenormalised re- 
sults can be produced automatically. 

For computations involving loop propagators for mas- 
sive fermions, we introduced the automatic generation 
of a mass counter term needed for the on-shcll renormal- 
isation of the massive particle. Here, we exploit the fact 
that such a counter term is strictly related to the mas- 
sive fermion self energy bubble diagrams (see Fig. [3]). As 




Fig. 3 Feynman diagram of a massive quark self energy in 
QCD. For this type of diagram GoSam automatically gener- 
ates UV-counterterms. 



described in Section |2"T21 the program GoSam analyzes 
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all generated diagrams. In that step also self-energy in- 
sertions of massive quarks are detected, where we make 
the replacement 



+ f + to) • 



'• + f + to) ■ gt*" 



[(q + r) 2 — to 2 ] q 2 [(q + r) 2 — to 2 ] q 2 



m 

T 



6q ■ r + 3(r 2 — to 2 ) 



3(4+1 H vV 2 



3to 2 



[(q + r) 2 — m 2 ] q 2 



(26) 



The symbol Ihv is one in the 't Hooft Veltman scheme 
and zero in DRED. 

Performing the integral, contracting the expression 
with the QCD vertices at both sides and multiplying 
the missing factor of (2tt)~ 1 we retrieve the expression 
for the mass counter-term, 



Sm 



(4tt) £ C f ( // 



27rr(l-e) 2 



5 - 1 



e 



HV 



(27) 



Furthermore, the renormalisation of a s leads to a 
term of the form 



\M\ 



2 

ct, ct s 



b- 



27rr(l - e) 



\M\ 



2 

tree 



e 



orr N,+N f , h 



3e 



q=N f +l 



+ %(1-1HV) 



(28) 



with /3 = (HCU - 4T R N f )/6, N f being the number of 
light quark flavours, Nf.h the number of heavy flavours, 
and b is the power of the coupling in the Born amplitude 
as defined in Eq. ([2~4"|). The last term of Eq. pro- 
vides the finite renormalisation needed to compensate 
the scheme dependence of a s , 



,DR 



1 + 



C A a 



MS 



6 2tt 



(29) 



A further contribution consists of the wave-function 
renormalisation of massive external quark lines. If we 
denote the set of external massive quark lines by Qh = 
{Qi(mi), . . . , Q p (m p )} we obtain 



\M\ 2 



a s (Any C t 



wf, 



27i-r(l-e) 2 



E 

Q(m)sQ h 



e 



lhv 



\M\ 



2 

tree 



(30) 



Finally, also the wave function of the gluon receives 
a contribution from the presence of heavy quarks in 



closed fcrmion loops. If N g is the number of external 
gluon lines, this contribution can be written as 



\M\L 



a s UttY Ar 2T R 
■N„ — 



27rr(l-e) 



3t 



Nf+N f:h 2 

E 

^ V TO 2 

q= Nf+l v 1 



1-MlLe. (31) 



At the level of the generated Fortran code the pres- 
ence of these contributions can be controlled by a set 
of variables defined in the module config.f90. The 
variable renormalisation can be set to 0, 1, or 2. 
If renormalisation=0, none of the counterterms ar 
present. If renormalisation=2 only |A^|^ t § m is in- 
cluded, which is the counterterm stemming from all 
terms of the type of Eq. (|2"T)) contributing to the am- 
plitude. 

In the case where renormalisation=l a more fine- 
grained control over the counterterms is possible. 

renorm_logs: if set to false, in all counterterms the 
generation of logarithms is disabled, i.e. factors of 
the form (») £ in eqs. ((27)) to (f3Tj) are replaced by 
one. 

renorm_beta: if set to false, the counterterm |.M| 2 t Qs 
is set to zero. 

renormjnqwf : if set to false, the counterterm |A^|^ f q 
is set to zero. 

r enormjnqs e : if set to false, the counterterm | Ai 1 2 t SnlQ 
is set to zero. 

renorm_decoupling if set to false, the counterterm 
|.M|^f g is set to zero. 

The default settings for renormalisation=l arc true 
for all the renorm options listed above. 



Finite Renormalisation of 75 in QCD In the 't Hooft 
Veltman scheme, a finite renormalisation term for 75 
is required beyond tree level. The relevant terms are 
generated only if fr5 is added in the input card to 
the list of extensions before code generation. Currently, 
the automatic generation of this finite contribution is 
not performed if model files different from the built-in 
model files are used. In agreement with |72j and |73j we 
replace the axial component at each vertex, 



7^75 -> ^axial (7^75 " 757^) , 



with 



Z; 



axial 



i 



2— -Cf ■ Ihv- 

Z7T 



(32) 



(33) 



Once it is generated, this contribution can be switched 
on and off at run-time through the variable renorm_gam- 
ma5, which is defined in the module config.f90. 
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Conversion between the Schemes In GoSam we have 
implemented two different schemes, the 'tHooft Velt- 
man scheme and dimensional reduction. By default, 
the former is used, while the latter can be activated 
by adding the extension dred. If a QCD computation 
has been done in dimensional reduction the result can 
be converted back to the 't Hooft Veltman scheme by 
adding a contribution for each external massless parton, 

I-m cdr ILoop = |A4 DR liioo P 

-Next 

-gl^ DR lt 2 r ccE7? R (34) 
1=1 

with 7™ = jf R = C F /2 and j® R = C A /6. This con- 
version can be switched on by setting convert_to_cdr 
to true in the module config.f90. At one-loop level, 
the 't Hooft Veltman scheme and conventional dimen- 
sional regularisation (CDR) are equivalent in the sense 
that 7^ tHV = for all partons. 

3 Requirements and Installation 

3.1 Requirements 

The program GoSam is designed to run in any mod- 
ern Linux/Unix environment; we expect that Python (> 
2.6), Java (> 1.5) and Make arc installed on the system. 
Furthermore, a Fortran 95 compiler is required in or- 
der to compile the generated code. Some Fortran 2003 
features are used if one wants to make use of the Les 
Houches interface [55] ■ We have tried all examples using 
gf ortran versions 4.1 and 4.5. 

On top of a standard Linux environment, the pro- 
grams FORM gB], version > 3.3 (newer than Aug. 11, 2010) 
and QGRAF [5] need to be installed on the system. Where- 
as spinney [49] and haggies [50] are part of GoSam 
and are not required to be installed separately, at least 
one of the libraries SAMURAI [40] and Golem95C [52] 
needs to be present at compile time of the generated 
code. Optionally, PJFRY [35] can be used on top of Go- 
lem95C. 



3.2 Download and Installation 

QGRAF The program can be downloaded as Fortran 
source code from 

http : //cf if . ist . utl . pt/~paulo/ qgraf . html . 

After unpacking the tar-ball, a single Fortran 77 file 
needs to be compiled. 



FORM The program is available at 

http : //www . nikhef . nl/~f orm/ 

both as a compiled binary for many platforms and as 
a tar-ball. The build process, if built from the source 
files, is controlled by Autotools. 

SAMURAI and Golem95C These libraries are available 
as tar-balls and from subversion repositories at 

http : //pro j ects . hepf orge . org/ samurai/ 

and 

http : //pro j ects . hepf orge . org/ golem/95/ 

respectively. For the user's convenience we have pre- 
pared a package containing Samurai and Golem95C 
together with the integral libraries OneLOop [71] . QCD- 
Loop |75| and FF |59j . The package gosam-contrib- 
l.O.tar.gz containing all these libraries is available 
for download from: 

http: //projects .hepf orge . org/gosam/ 

GoSam The user can download the code either as a 
tar-ball or from the subversion repository at 

http : //pro j ects . hepf orge . org/gosam/ . 

The build process and installation of GoSam is con- 
trolled by Python Distutils, while the build process 
for the libraries Samurai and Golem95C is controlled 
by Autotools. 

Therefore the installation proceeds in two steps: 

1. For all components which use Autotools, the fol- 
lowing sequence of commands installs them under 
the user defined directory MYPATH. 

./configure --pref ix=MYPATH 

make FC=gfortran F77=gfortran 

make install # or sudo make install 

If the configure script is not present, the user needs 
to run sh . / autogen . sh first . 

2. For GoSam which is built using Distutils, the 
user needs to run 

python setup. py install \ 
--prefix MYPATH 

If MYPATH is different from the system default (e.g. 
/usr/bin), the environment variables PATH, LD_LIB- 
RARY_PATH and PYTH0NPATH might have to be set 
accordingly. For more details we direct the user to 
the GoSam reference manual and to the documen- 
tation of the beforcmcntioncd programs. 
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4 Using GoSam 

4.1 Setting up a simple Process 

GoSam is a very flexible program and comes with a 
wide range of configuration options. Not all of these 
options are relevant for simple processes and often the 
user can leave most of the settings at their default val- 
ues. In order to generate the code for a process, one 
needs to prepare an input file, which will be called pro- ^ 

o 

cess card in the following, which contains 

— process specific information, such as a list of initial 
and final state particles, their helicities (optional) 
and the order of the coupling constants; 

— scheme specific information and approximations, such 
as the regularisation and renormalisation schemes, 
the underlying model, masses and widths which are 
set to zero, the selection of subsets of diagrams; the 
latter might be process dependent; 

— system specific information, such as paths to pro- 
grams and libraries or compiler options; 

— optional information for optimisations which control 
the code generation. 

In the following we explain how to set up the required 
files for the process qq — > gZ° — > ge~e + . The example 
computes the QCD corrections for the uu initial state, 
where m e = and Nf = 5 massless quarks are as- 
sumed. For our example, we follow an approach where 
we keep the different types of information in separate 
files - process. rc, scheme. rc and system. rc - and 
use GoSam to produce a process card for this process 
based on these files. This is not required — one could 
also produce and edit the process card directly — it is 
however more convenient to store system specific infor- 
mation into a separate, re-usable file, and it makes the 
code generation more transparent. 

Process specific information The following listing con- 
tains the information which is specific to the process. 
The syntax of process cards requires that no blank char- 
acter is left between the equals sign and the property 
name. Commentary can be added to any line, marked 
by the '#' character. Line continuation is achieved using 
a backslash at the end of a line! 

Listing 1 File 'process. rc' 

1 process_path=qqgz 

2 in=u , u~ 

3 out=g , e- , e+ 

4 helicities=+-+-+ ,+-- 
•5 order = QCD , 1 , 3 



The first line defines the (relative) path to the direc- 
tory where the process files will be generated. GoSam 
expects that this directory has already been created. 
Lines 2 and 3 define the initial and final state of the 
process in terms of field names, which arc defined in 
the model file. Besides the field names one can also use 
PDG codes [751177] instead. Hence, the following lines 
would be equivalent to lines 2 and 3 in Listing [TJ 




Line 4 describes the helicity amplitudes which should 
be generated. If no helicities are specified, the program 
defaults to the generation of all possible helicity con- 
figurations, some of which may turn out to be zero. 
The different helicity amplitudes are separated by com- 
mas; within one helicity amplitude there is one charac- 
ter (usually '+', '-' and '0') per external particle from 
the left to the right. In the above example for the reac- 
tion 

u(fci,Ai)u(fc 2 , A 2 ) -> .g(fc 3 , A 3 )e~(fc 4 , A 4 )e + (fc 5 ,A 5 ) 
we have the following assignments: 



Helicity 


Ai A 2 A3 A4 A5 





+ - + - + 


1 


+ - - - + 


2 


+ + - + 


3 


+ - - + 



With the above value for helicities we generate all 
non-vanishing helicities for the partons but keep the 
lepton helicities fixed. In more complicated examples 
this way of listing all helicities explicitly can be very 
tedious. Therefore, we introduced the option to gener- 
ate sets of helicities using square brackets. For example, 
if the gluon helicity is replaced by [+-] , the bracket is 
expanded automatically to take the values + ,-. 

helicities = +- [+-]-+ , -+[+-]- + 

A further syntactical reduction can be achieved for the 
quarks. The current expansion of a square bracket and 
its opposite value can be assigned to a pair of variables 
as in [xy=+-] . If the bracket expands to '+' then x is 
assigned '+' and y is assigned the opposite sign, i.e. 
If the bracket expands to '-' the assignments are x=- 
and y=+. Hence, the helicity states of a massless quark 
anti-quark pair are generated by [qQ=+-] Q, and the se- 
lection of helicities in our example can be abbreviated 
to 



-+ , -++-+ , -+--+ 

4 helicities= [qQ=+-] Q [+-] -+ 



5 The line numbers are just for reference and should not be 
included in the actual files. 



which is equivalent to the version of this line in List- 
ing [TJ 



12 



Finally, the order (power) of the coupling constants 
has to be specified. Line 5 contains a keyword for the 
type of coupling (QCD or QED), the order of this coupling 
constant in the unsquared tree level amplitude (in our 
example: 1) and the order of the coupling constant in 
the unsquared one- loop amplitude (in our example: 3). 
One can also use GoSam to generate the tree level only, 
by giving only the power of the tree level amplitude: 

5 order=QCD,l 

Conversely, GoSam will generate the virtual amplitude 
squared for processes where no tree level is present if 
the tree level order is replaced by the keyword NONE. 

5 order = QCD , NONE ,3 

Up to now, the file would generate all 8 tree level 
and 180 one-loop diagrams contributing to the process 
uu — > ge~e + , regardless of the intermediate states. 
Nevertheless, what we intended to generate were only 
those diagrams where the electron pair comes from the 
decay of a Z — >• e~e + . GoSam offers two ways of achiev- 
ing this diagram selection, cither by passing a condition 
to QGRAF or by applying a filter written in Python. The 
first option would be specified by the option qgraf . 
verbatim, which copies the argument of the option to ^ 

o 

the QGRAF input file in verbatim. The following filter 
demands the appearance of exactly one Z-propagator, ^ 



1 



leaving us with 2 tree-level and 45 one-loop diagrams: 
6 qgraf . verbat im= true= iprop [Z , 1 , 1] ; 



G 
7 

8 

The alternative solution is the application of a Python 
filter using the options filter, lo for tree level and ^ 
filter. nlo for one-loop diagrams. The current exam 
pie requires the two lines 



6 filter. lo= IPROP ( [Z ])== 1 

7 filter. nlo= IPROP ( [Z ])== 1 
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Scheme specAfic information For our example we put 
all scheme specific definitions in the file scheme. rc. It 
contains the choice of a suitable regularisation scheme 
and fixes what types of UV counterterms are included 
in the final result. 

Listing 2 File 'scheme. rc' 

1 extensions=dred 

2 qgraf . opt ions= onshell 

3 zero=mU , mD , mC , mS , mB , me , wT 

4 one=gs 

In Listing [2j line 1 selects dimensional reduction as 
a regularisation scheme. If dred is not specified in the 
list of extensions, GoSam works in the 't Hooft Veltman 
scheme by default. Line 2 removes all on-shell bubbles 




on external legs. This is, on the one hand, required to 
be consistent with the renormalisation scheme. On the 
other hand, those diagrams would lead to zero denom- 
inators at the algebraic level. In line 3 all light quark 
masses, the mass of the electron and the width of the 
top quark are set to zero. Further, as a convention 
rather than a scheme, the strong coupling g s is set to 
one in line 4, which means that g s will not occur in the 
algebraic expressions, assuming that the user will mul- 
tiply the final result by his desired value for the strong 
coupling. If the option one=gs is not used, the default 
value contained in the file common/model . f 90 will be 
used. This default value of course can be changed by 
the user. 

System specific information In order to adapt the code 
generation to the system environment, GoSam needs to 
find a way of determining all relevant paths and options 
for the programs and libraries used during generation, 
compilation and linking of the code. Those settings are 
fixed in the file system. rc in our example^ 

Listing 3 File 'system. rc' 

system . ext ens ions= samurai , golem95 
samurai . fcflags=\ 

-1$ {PREFIX}/ include/ samurai 
samurai . ldflags=\ 

-L${PREFIX}/lib -lsamurai 
samurai . version=2 . 1 . 1 
golem95 . fcflags=\ 

- 1$ {PREFIX}/ include/ golem95 
golem95. ldflags=\ 

-L${PREFIX}/lib -lgolem95 
f orm . bin = ${ PREFIX }/ bin/ 1 form 
qgraf . bin=${PREFIX}/bin/qgraf 
fc . bin=gfortran 

Generating the Code After having prepared the input 
files correctly we need to collect the information dis- 
tributed over the three files process. rc, scheme. rc 
and system. rc in one input file, which we will call 
gosam. in here. The corresponding command is: 

gosam.py --template gosam. in \ 
--merge process. rc \ 

--merge scheme. rc --merge system. rc 

The generated file can be processed with gosam.py di- 
rectly but requires the process directory to be present. 

mkdir qqgz 
gosam.py gosam. in 
cd qqgz 

6 In this example we assume that the user has defined an 
environment variable PREFIX. 
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All further steps are controlled by the generated make 
files; in order to generate and compile all files relevant 
for the matrix element one needs to invoke 

make compile 

The generated code can be tested with the program 
matrix/test . f 90. The following sequence of commands 
will compile and run the example program. 

cd matrix 
make test . exe 
. /test . exe 



The last lines of the program output should look as 
follow^ 



# 






L0 : 


. 3450350717601E-06 


# 


NLO , 


f inite 


part 


-10 . 77604823456547 


# 


NLO , 


single 


pole 


-19 . 98478948141949 


# 


NLO , 


double 


pole 


-5 . 666666665861926 


# 


IR , 


single 


pole 


-19 . 98478948439310 


# 


IR , 


double 


pole 


-5 . 666666666666666 



The printed numbers are, in this order, ao, cq/oq, C-x/oq, 
c -2/ao and the pole parts calculated from the infrared 
insertion operator [751175] , 

One can generate a visual representation of all gen- 
erated diagrams using the command 

make doc 

which generates the file doc/process .ps using a Py- 
thon implementation of the algorithm described in [60) 
and the M^X package AXODRAW gT]. 

4- 1.1 Further Options 

GoSam provides a range of options which influence 
the code generation, the compilation and the numer- 
ical evaluation of the amplitude. Giving an exhaustive 
list of all options would be far beyond the scope of this 
article and the interested user is referred to the refer- 
ence manual. Nonetheless, we would like to point out 
some of GoSam's capabilities by presenting the corre- 
sponding options. 

Generating the R2 Term When setting up a process 
the user can specify if and how the R2 term of the 
amplitude should be generated by setting the variable 
r2 in the setup file. 

r2=explicit 

7 The actual numbers depend on the random number genera- 
tor of the system because the phase space point is generated 
randomly; however, the pole parts should agree between the 
matrix element and the infrared insertion operator given that 
the matrix element is fully renormalised. 



Possible options for r2 are implicit, which is the de- 
fault, explicit, off and only. The keyword implicit 
means that the R2 term is generated along with the 
four-dimensional numerator as a function in terms of q, 
/i 2 and e and is reduced at runtime by sampling differ- 
ent values for /i 2 . This is the slowest but also the most 
general option. Using the keyword explicit carries out 
the reduction of terms containing fj? or e during code 
generation (see Appendix [Bj . The keyword off puts 
the i?2 term to zero which is useful if the user wants to 
provide his own calculation for these terms. Conversely, 
using r2=only discards everything but the R2 term (re- 
ducing it as in the case explicit) and puts GoSam in 
the position of providing R2 terms for external codes 
which work entirely in four dimensions. 

Diagram Selection GoSam offers a two-fold way of se- 
lecting and discarding diagrams. One can either influ- 
ence the way QGRAF generates diagrams or apply fil- 
ters to the diagrams after they have been generated by 
QGRAF or combine the two methods. Let us assume that 
in the above example we want to remove the third gen- 
eration of quarks completely. Hence, all closed quark 
loops would be massless and therefore the second gen- 
eration is just an exact copy of the first one. We can 
therefore restrict the generation of closed quark loops 
to up and down quarks. GoSam has a filter precisely for 
this purpose, which takes the field names of the flavours 
to be generated as arguments. 

filter .nlo=NFGEN(U,D) 

This filter can be combined with the already existing 
filter selecting only diagrams containing a Z-propagator 
using the AND function: 

filter .nlo=AND( NFGEN (U , D ) , \ 

IPRDP([Z]) == 1 ) 

A further feature of the code generated by GoSam 
is the possibility of selecting diagrams at runtime. For 
example, we would like to distinguish at runtime three 
different gauge invariant sets of diagrams at one-loop 
level: 

1. diagrams with a closed quark loop where the Z is 
attached to the loop; 

2. diagrams with a closed quark loop where the Z is 
emitted from the external quark line; 

3. diagrams without a closed quark loop. 

In order to provide the code for a diagram selection at 
runtime one simply replaces the above filter by a list of 
filters as follows 

filter . nlo= [\ 

AND ( NFGEN (U , D ) , IPROP ( [Z] ) == 1, \ 
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NF, LODPVERTICES ( [Z] , _ , _) ==1), \ 
AND ( NFGEN (U , D) , IPROP ( [Z] ) == 1, \ 

NF, LODPVERTICES ( [Z] ,_, _) ==0), \ 
AND ( NFGEN (U , D) , IPROP ( [Z] ) == 1, \ 

NOT ( NF ) ) ] 

The two new filters in use are NF which selects closed 
quark loops only and LODPVERTICES which counts the 
number of vertices attached to the loop with the given 
sets of fields running through the vertex (where _ re- 
places any held). In the Fortran files one can access the 
diagram selection through the routine update Jlags. 
The three selection criteria are stored in a derived data 
type virt Jlags which has fields eval_0, . . . , eval_2, 
in general ranging from zero to the length of the list 
given in filter. nlo. 

use groups 

type (virt_f lags) :: flags 

f lags%eval_0= . true . ! first group only 

flags7 eval_l = . false . 

flags7 eval_2=. false . 

call updat e_f lags ( f lags ) 

Additional Extensions Some of GoSam's functionality 
is available through the extensions variable. On top 
of the already presented options for selecting a regular- 
isation scheme (by adding the option dred) or for acti- 
vating interfaces to several different reduction libraries 
(samurai, golem95, pjfry) the user can also add the 
following options: 

fr5 adds and activates the relevant code for the com- 
putation of the finite renormalisation of 75 required 
in the 't Hooft Veltman scheme as described in Eq. (1321) . 

powhegbox generates routines for the computation of 
the color and spin correlated Born matrix elements 
as required by POWHEG [34] . 

autotools uses make files which use Autoconf and Au- 
tomake for compilation of the matrix element. 

gauge check replaces the polarisation vectors of exter- 
nal vector fields by 

e"(fc i )^e"(fc i ) + 2i^ (35) 

where the variable Zj defaults to zero and is ac- 
cessible in the Fortran code through the symbol 
gauge iz. 

4.2 Interfacing the code 

The matrix element code generated by GoSam pro- 
vides several routines to transparently access partial or 
full results of the amplitude calculation. Here, we only 



present a minimal set of routines which can be used 
to obtain the set of coefficients [a$, Co, c_i, c_2] for a 
given scale and a given set of external momenta. The 
routines, which can be accessed through the modules 
matrijo arc defined as follows: 

initgolem This subroutine must be called once before 
the first matrix element evaluation. It initializes all 
dependent model parameters and calls the initiali- 
sation routines of the reduction libraries. 

interface 

subroutine initgolem( init_libs) 

use config, only: ki 

logical, optional, & 
& intent(in) :: init_libs 

end subroutine 
end interface 

The optional argument init_libs can usually be 
omitted. It should be used only when several ini- 
tialisation calls become necessary, but the reduction 
libraries and loop libraries should be initialized only 
once. All model parameters are accessible as global 
variables in the module model and should be modi- 
fied (if at all) before calling initgolem. 
Samplitude This subroutine starts the actual calcula- 
tion of the amplitude for a given phase space point. 



interface 




subroutine Samplitude 


t 


& (vecs , scale2 , 


amp , ok , h) 


use config, only: ki 




use kinematics, only: 


num_legs 


real(ki) , dimension(num_legs ,4) ,& 


& intent (in) : 


: vecs 


double precision, & 




& intent (in) : 


: scale2 


double precision, & 




& intent ( out ) 


: : amp 


logical, optional, & 




& intent ( out ) 


: : ok 


integer, optional, & 




& intent (in) : 


: h 


end subroutine 




end interface 





The first mandatory arguments of this routine are 
the external momenta vecs, where vecs (i , : ) con- 
tains the momentum of the i-th particle as a vec- 
tor [Ei,pf ,p\ ,pf], and we use in-out kinematics, i.e. 
Pi + Pi — > P3 + • • ■ + Pn- Maximal numerical sta- 
bility is achieved if the beam axis is chosen along 

8 If a process name was given all modules are prefixed by the 
name, e.g. if process_name=pr01, the module matrix would be 
renamed into prOljnatrix. 
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the z-axis. The second argument, scale2 = fj^, is 
the square of the rcnormalisation scale. As a third 
argument the routine expects a vector which ac- 
cepts the result in the format [ao, Co, c_i, C-2] with 
the cocfficents being defined in Eqs. (|24l) and (|25j) . 
The optional argument ok may be used in order to 
report the outcome of the reconstruction tests in 
samurai if no rescue method has been chosen (see 
Section 12.4. 2\i . The last argument allows one to se- 
lect a single helicity subamplitude; the index h runs 
from zero to the number of helicities minus one. The 
labeling of the helicities is documented for each pro- 
cess in the file doc/process . ps. 
exitgolem This routine should be called once after the 
last amplitude evaluation in the program. It closes 
all open log files and gracefully terminates the re- 
duction and loop libraries. 

interface 

subroutine exitgolem ( exit _1 ibs ) 

use config, only: ki 

logical, optional, 
& intent(in) :: exit_libs 

end subroutine 
end interface 

The optional argument exit_libs should only be 
set if multiple calls to this routine (e.g. for different 
matrix elements) are necessary and the dependent 
libraries should be terminated only once. 

A small program which computes the amplitude for a 
set of phase space points is automatically generated 
with the amplitude code in the file test. f 90 in the 
subdirectory matrix. The script config. sh in the pro- 
cess directory returns suitable compilation and linking 
options for the generated matrix element code. 



4.3 Using the BLHA Interface 

The so-called Binoth Les Houches Accord (BLHA) [56] 
defines an interface for a standardized communication 
between one-loop programs (OLP) and Monte Carlo 
(MC) tools. The communication between the two sides 
is split into two main phases: an initialisation phase and 
a runtime phase. During initialisation the two programs 
establish an agreement by exchanging a set of files and 
typically initiate the code generation. The OLP runtime 
code is then linked to the MC program and, during the 
runtime phase, called through a well-defined set of rou- 
tines providing NLO results for the phase space points 
generated by the MC. According to this standard, it 
is the responsibility of the MC program to provide re- 
sults for the Born matrix element, for the real emission 



wrile order file 



ear] eonlrael Trie 



read order file 



wrile eonlrael file 



eall OLP^Ev; 
Subprooess 



Code Generation and Lnikine 



initialization 
phase 



initialize OLP 
compute result 
return result 



Fig. 4 Schematic overview over the interaction between 
Monte Carlo tool and one-loop program in the Binoth Les 
Houches Accord. 



and for a suitable set of infrared subtraction terms. A 
schematic overview on this procedure is shown in Fig. [4] 
GoSam can act as an OLP in the framework of the 
BLHA. In the simplest case, the MC writes an order 
file — in this example it is called olp_order . In — and 
invokes the script gosam.py as follows: 

gosam.py --olp olp_order.lh 

Further, GoSam specific options can be passed either 
in a file or directly at the command line. One can, for 
example, use autotools for the compilation by modify- 
ing the above line as follows. 

gosam.py --olp olp_order . In \ 
extensions= autotools 

The contract file is given the extension .olc by de- 
fault and would be olp_order . olc in this example. Al- 
ternatively, the name can be altered using the -o op- 
tion. 

If successful, the invocation of gosam.py generates 
a set of files which can be compiled as before with a 
generated make file. The BLHA routines are defined 



in the Fortran module olpjnodule but can also be 
accessed from C programrjj. The routines 0LP_Start 
and OLP_EvalSubProcess are defined exactly as in the 
BLHA proposal [56] . For convenience, we extended the 
interface by the functions 0LP_Finalize () , which ter- 
minates all reduction libraries, and 0LP_0ption(char* , 
int*) , which can be used to pass non-standard options 
at runtime. For example, a valid call in C to adjust the 
Higgs mass would be 



int 


ierr ; 






OLP 


.Options ("mH = 146 . 78" , 


& ierr ) ; 





A value of one in ierr indicates that the setting was 
successful. A value of zero indicates an error. 



3 A header file is provided in olp.h. 
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4.4 Using External Model Files 

With a few modifications in the process description 
files, GoSam can immediately make use of model files 
generated by either FeynRules [5U] in the UFO format [53] 
or by LanHEP [55] ■ In both cases, the following limita- 
tions and differences with respect to the default model 
files, sm and smdiag, apply: 

— As usual, particles can be specified by their PDG 
code. The field names, as used by QGRAF, are part?' 
and antii for the particles with the PDG code i and 
—i respectively. For example, the W + and the W~ 
boson would be called part24 and anti24. 

— All model parameters are prefixed by the letters mdl 
in order to avoid name clashes with existing variable 
names in the matrix element code. 

— The variable model . options and the extension f r5 
are not guaranteed to work with models other than 
the built-in models. 

Importing models in the UFO format A model descrip- 
tion in the UFO format consists of a Python package 
stored in a directory. In order to import the model into 
GoSam one needs to set the model variable in the input 
card to specify the keyword FeynRules in front of the 
directory name, where we assume that the model de- 
scription is in the directory $H0ME/models/MSSM_UF0. 

model=FeynRules , $H0ME/models /MSSM_UF0 

Importing models in the LanHEP format LanHEP mo- 
del descriptions consist of a set of plain text files in 
the same directory with a common numbering (such as 
func4.mdl, lgrng4.mdl, prtcls4.mdl, vars4.mdl). A 
LanHEP model can be loaded by specifying the path 
and the common number in the model variable. As- 
suming the files are situated in the directory $H0ME/ 
models/MSSM_LHEP one would set the variable as fol- 
lows. 

model=$HOME/models/MSSM_LHEP ,4 

Details about the allowed names for the table columns 
are described in the GoSam reference manual. Precom- 
piled MSSMJJFO and MSSM_LHEP files can also be found 
in the subdirectory examples/model. 

5 Sample Calculations and Benchmarks 

The codes produced by GoSam have been tested on 
several processes. In this section we describe some ex- 
amples of applications. Additional results, whose corre- 
sponding code is also included in the official distribution 
of the program, will be reported in Appendix B. 



5.1 pp ->W~+j with SHERPA 

In Section |4~51 the BLHA interface of GoSam was pre- 
sented. This interface allows one to link the program 
to a Monte Carlo event generator, which is, in gen- 
eral, responsible for supplying the missing ingredients 
for a complete NLO calculation of a physical cross sec- 
tion. Among the different general purpose Monte Carlo 
event generators, SHERPA [81] is one of those which 
offers these tools: computing the LO cross section, the 
real corrections with both the subtraction terms and 
the corresponding integrated counterparts [5211551154"] . 
Furthermore, SHERPA offers the possibility to match 
a NLO calculation with a parton shower [551155] . Using 
the BLHA interface, we linked GoSam with SHERPA 
to compute the physical cross section for W~ + 1 jet at 
NLO. 

The first steps to perform this linking is to write a 
SHERPA input card for the desired process. Instruc- 
tions and many examples on how to write this can be 
found in the on-line manual [87] . Running the code for 
the first time will produce an order file OLE_order.lh 
which contains all the necessary information for Go- 
Sam, to produce the desired code for the loop part of 
the process. This includes a list of all partonic sub- 
processes needed. In parallel to the production of the 
needed SHERPA libraries with the provided script, one 
can at this point run the gosam.py command with the 
flag — olp and the correct path to the order file as 
explained in Section |4~51 Further options may be spec- 
ified. Among them it is useful to have a second, Go- 
SAM-specific, input card with all the important GoSam 
options. Since, at the end, SHERPA needs to be linked 
to a dynamic library, it is convenient to run GoSam 
with the autotools extension, which allows the direct 
creation of both static and dynamic libraries, together 
with the test routine test. The gosam.py script cre- 
ates all the files needed for interfacing GoSam with the 
Monte Carlo event generator together with the code for 
the one-loop computation of all needed subprocesses, 
and a makefile to run them. The different parton-level 
subprocesses are contained in different subdirectories. 
At this point the user simply has to run the makefile to 
generate and compile the code. Once the one-loop part 
of the code is ready, the produced shared library must 
be added to the list of needed libraries in the SHERPA 
input card as follows. 

SHERP A_LD ADD = LHOLE golem_olp; 

With this operation the generation of the code is com- 
pleted. The evaluation of the process and the physical 
analysis can then be performed at the user's discretion 
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following the advice given in the SHERPA on-line doc- 
umentation [57]. 

We tested the BLHA interface by computing W~ + 1 
jet and producing distributions for several typical ob- 



parameters 



servables. In Figs. 5(a) and 5(b) the inclusive transverse 



momentum and rapidity of the jets is shown. These dis- 
tributions were compared with similar ones produced 
using the program MCFM [2"Tll28] . and perfect agree- 
ment was found. 



5.2 pp ->■ W ± + j, EW Corrections 

As a first example of an clectroweak calculation, we 
computed the virtual one-loop corrections to ud — ¥ Wg. 
A complete analytical calculation for this process was 
presented in Ref. [55] , 

parameters 



M z 


91.1876 


M w 


80.419 


cos W 


0.88156596117995394232 


fj, 


M w 



For the kinematic point given in Tab. [2] and the above 
parameters we obtain the following result: 







result ud — 


>Wg 


co /ao unren. 
c_i/ao unren. 
c-2/ao unren. 




2.812364835883295 
-94.52525523327047 
17.84240236996827 
-0.5555555555555560 






renormali 


zed 




GoSam 




Eqs. (67,70) of Ref. [HS] 


c_i/a 
c-2/ao 


4.743825167813529 
-0.5555555555555560 


4.7438251678146885 
-0.5555555555555555 



The poles have been renormalizcd using Eqs. (49)- 
(64) in Sections 3.3 and 3.4 of [EE]- Our result is agree- 
ment with Eqs. (67), (70) of Ref. [EE] and with Ref. [EH] 
for the infrared divergences that remain after renormal- 
isation. 



5.3 77 — > 77 

The process 77 — > 77 in the Standard Model first arises 
at the one-loop order, and proceeds through a closed 
loop of fermions and W bosons. Of the 16 helicity am- 
plitudes contributing to it, only three are independent 
and their analytic expressions can be found in [5D] . The 
pure QED contribution, involving a fermion loop, is 
contained in samurai-1 . |40] and will not be repeated 
here. Instead, we show the results of the W-loop con- 
tribution to the independent helicity amplitudes, as an 
example of EW corrections that can be handled with 
GoSam. 



y/s 1000 


/' 




M w 80.376 


e 


1 



With the above parameters and the kinematics of Tab. [3] 
we obtain the following results. 



result 77 — > 77 (EW) 




GoSAM(dred) 


Refs.jM] 




12.02541904626610 


12.025419045962 


\M++-+\ 


7.380406043429961 


7.3804060437434 


\M + + —\ 


982.7804939723322 


982.78049397093 



5.4 pp -> X1X1 m the MSSM 

As an example for the usage of GoSam with a model file 
different from the Standard Model we calculated the 
QCD corrections to neutralino pair production in the 
MSSM. The model file has been imported via the in- 
terface UFO (Universal FeynRules Output) [S3] which 
facilitates the import of Feynman rules generated by 
FeynRules [80] to programs generating one-loop am- 
plitudes. To import such files within the GoSam setup, 
all the user has to do is to give the path to the corre- 
sponding model file in the input card. 

For this example, we combined the one-loop am- 
plitude with the real radiation corrections to obtain 
results for differential cross sections. A calculation of 
neutralino pair production for the LHC presenting to- 
tal cross sections at NLO is given in [HI] . 

For the infrared subtraction terms the program Mad- 
Dipole [92 , 933 is used, the real emission part is calcu- 
lated using MadGraph/MadEvent [M]- The virtual ma- 
trix element is renormalized in the MS scheme, while 
massive particles are treated in the on-shell scheme. 
The renormalisation terms specific to the massive MSSM 
particles have been added manually. 

In Fig. [5] we show the differential cross section for 
the fn x o x o invariant mass, where we employed a jet veto 
to suppress large contributions from the channel qg — > 
X1X1Q which opens up at order a 2 a s , but for large pi^ 
belongs to the distinct process of neutralino pair plus 
one hard jet production at leading order. We used Nf = 
5 massless quark flavours and the MSTW08 [HS] parton 
distribution functions. For the SUSY parameters we use 
the modified benchmarks point SPSlamod suggested in 
[9"B] , and we use yfs = 7 TeV. 

For reference, we also give the result for the unrenor- 
maliscd amplitude at one specific phase space point for 
uu X1X1 i n t nc DRED scheme, using the following 
parameters and momenta: 
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Fig. 5 NLO calculation of W +1 jet production at LHC using GoSam interfaced with SHERPA via the BHLA interface. The 
comparison to MCFM is also shown. 





E 


Px 


Py 


Pz 


11 


500 








500 


d 


500 








500 


W 


503.23360778049988 


110.20691318538486 


441.95397288433196 


-198.26237811718670 


a 


496.76639221950012 


-110.20691318538488 


-441.95397288433202 


198.26237811718664 



Table 2 Kinematic point used in pp — > + j, EW. 





E 


Pre 


Py 


Pz 


7 


500 








500 


7 


500 








-500 


7 


500 


436.6186300198938284 


-59.1784256571505765 


236.3516148798047425 


7 


500 


-436.6186300198938284 


59.1784256571505765 


-236.3516148798047425 



Table 3 Kinematic point used in 77 — > 77. 



parameters 



M z 


91.1876 


r z 





M w 


79.829013 


sin 2 9 W 


1 - M^/M 2 


fj, 


M z 


Nf 


1 


9s 


1 


a^ 1 


127.934 




96.6880686 


M S 


607.713704 


Mu L 


561.119014 


Mu R 


549.259265 


M ho 


110.899057 


M Ho 


399.960116 



All widths have been set to zero; for further set- 
tings we refer to the model parameter files contained in 
the subdirectory examples/model/MSSM_UF0. We have 
checked that the pole terms of the renormalised ampli- 
tude cancel with the infrared poles from MadDipolc. 
For the phase space point given in Tab. [4] we obtain the 
following numbers. 

GoSam result uu — > X1X1 
a 0.8680577964243597- 10" 3 

c /a -31.9136615197871 
c_i/a 13.4374663711899 
c_ 2 /ao 2.6666666666667 



5.5 e+e -> e + e _ 7 in QED 

As an example of a QED calculation, we compared the 
virtual QED corrections for the process e + e~ — > e + e~"f 
with the results provided in |97j . The results compared 
in the table are the bare unrenormalised amplitudes 
in the 't Hooft Veltman scheme. No counterterms or 
subtraction terms have been added to the result. 



parameters 





1.0 


a 


7.2973525376 


10~ 3 


/' 




m e 


0.51099891 


10~ 3 



Using the parameters given above and the kinemat- 
ics of Tab. [5] we obtain the following results. 





result e+e — > e 


^e 7 




GoSam 


Ref. [97] 


a 


0.7586101468103622 


0.7586101468103619 


c /a 


0.5005827938274887 


0.5005828268263969 


c-x/ao 


0.0474506407008029 


0.0474506427003504 


c-2/aa 
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E 


Px 


Py 


Pz 


u 
u 

a 
Xi 

Xi 


1000 
1000 
1000 
1000 





42.3752677206678996 
-42.3752677206678996 





115.0009952646289548 
-115.0009952646289548 


1000 
-1000 

987.7401101322898285 
-987.7401101322898285 



Table 4 Kinematic point used in pp — > x?Xi m the MSSM. 



(in) 
(in) 
(out) 
(out) 



0.5 
0.5 

0.1780937847558600 
0.3563944406457374 
0.4655117745984024 



Px 




0.1279164180985903 
0.02860530642319879 
0.1565217245217891 



Py 

0.4999997388800458 
-0.4999997388800458 
0.05006809884093004 
0.1832142729949070 
0.1331461741539769 



Pz 




0.1133477415216646 
0.3043534176228102 
0.4177011591444748 



Table 5 Kinematic point used in e^ 



x-io- 6 




1600 



m(N11,N12) [GeV] 

Fig. 6 Comparison of the NLO and LO m x a x o distributions 

for the process pp —>■ x?X? with a jet veto on jets with p^f* > 
20 GeV and r? < 4.5. The band gives the dependence of the 
result on = up = ur between uq/2 and 2ug. We choose 
fj,o = Mz- The black line gives the bin error for the value at 
the central scale. 

5.6 pp ttH 

This process has been compared with the results given 
in [35] . The partonic subprocesses uu — > ttH and gg — > 
ttH where computed both in the 't Hooft Veltman scheme 
and in dimensional reduction and the fully renormalised 
results were successfully compared as an internal consis- 
tency check. Apart from wave function renormalisation 
and mass countertcrms, Yukawa coupling renormalisa- 
tion is also needed here. Yukawa coupling countertcrms 
are in this case equal to the wave function countertcrms. 
The Yukawa top mass is set equal to its pole mass. 



parameters 

Vs 500.0 Nf 5 

\i m t N fyh 1 

m t 172.6 a s 0.1076395107858145 

m H 130 v 246.21835258713082 



The kinematics used to obtain the results below is 
given in Tab. H3 The results are given in the 't Hooft 
Veltman scheme, and are fully renormalised. 





result uu — > 


ttH 




GoSam 


Ref. [39] 


a ■ 10 4 


2.200490364806190 


2.2004904613782828 


c /a 


-15.29615178164782 


-15.29615211731521 


c-i/a 


-1.640361500121837 


-1.640361536072381 


c-2/ao 


-2.666666666666666 


-2.666666725182165 




result gg — > ttH 




GoSam 


Ref. EH] 


a ■ 10 5 


6.127399805961155 


6.127400074872043 


eo/ao 


9.006680638719660 


9.006680836410272 


e_i/a 


2.986347664537282 


2.9863477301662056 


c-2/ao 


-6.000000000000004 


-6.000000131659877 



On an Intel Core i7 950 at 3 GHz the evaluation of 
a single phase space point took 44 ms in the uu channel 
and 223 ms in the gg channel. The code was compiled 
with gf ortran without optimisations. 



5.7 gg ttZ 

This amplitude, fully renormalised, has been compared 
with the results given in |37j . 



parameters 



9s 


1 


Gf 


0.0000116639 


/' 


m t 


N f 


5 


«'/ 


170.9 


M w 


80.45 


M z 


91.18 







The kinematics used to obtain the results below is given 
in Tab. El 



20 



u/g 

H 

t 
t 



250.0 
250.0 

136.35582793693018 
181.47665951104506 
182.16751255202476 



Px 

0.0 
0.0 

15.133871809486299 
20.889486679044587 
-36.023358488530903 



Py 

0.0 
0.0 

27.986733991031045 
-50.105625289561424 
22.118891298530357 



Pz 

250.0 
-250.0 

26.088703626953386 
14.002628607367491 
-40.091332234320859 



Table 6 Kinematic point used in pp — > ttH. 





E 


Px 


Py 


Pz 


g 


7000.0 


0.0 


0.0 


7000.0 


g 


7000.0 


0.0 


0.0 


-7000.0 


t 


6270.1855170414337 


-4977.7694025303863 


806.93726196887712 


3725.2619580634337 


i 


6925.5258180925930 


5306.3374282745517 


-1281.8763412410237 


-4258.3185872039012 


Z 


804.28866486597315 


-328.56802574416463 


474.93907927214622 


533.05662914046729 



Table 7 Kinematic point used for gg — > ttZ. 





result gg -» 


tiZ 




GoSam 


Ref. [37] 


a ■ 10 6 


0.1531395190212139 


0.1531395190212831 


Co /ao 


-204.9208290898557 


-204.920829867328 


c-i/ao 


50.62939646427283 


50.6293965717156 


c-2/ao 


-5.999999999999997 


-6.00000000000003 



The evaluation of a single phase space point took 1433 
ms on a 2 GHz processor. The code was compiled with 
gfortran -02. 





result uu — > 


bbbb 




GoSam 


Ref. [38] 


a ■ 10 9 


5.753293428094349 


5.753293428094391 


c /a 


-22.19223384585620 


-22.19223384564902 


c_i/a 


-20.89828996870689 


-20.89828996857439 


c-2/ao 


-8.000000000000199 


-8.000000000000037 



On an Intel Xeon E7340 the running time for the 
calculation of a single phase space point was 19.6 s for 
the gluon initiated channel and 440 ms for the quark 
channel. 



5.8 pp -> bbbb + X 

A detailed discussion of this process can be found in 
[98,99]. In this section we focus on the parts that are 
relevant in the context of the virtual corrections. In 
particular we compared our result to the one given in 
|38j , which is the fully renormalised amplitude including 
the mass counterterms for the top-quark contribution. 



parameters 





500 


Nf 


5 


/' 






1 




174 


m b 





r, 





9s 


1 



The results below are obtained for the phase space point 
of Tab. [5] using the above parameters. 



result gg — > bbbb 





GoSam 


Ref. [38j 


a ■ 10 6 


1.022839601391936 


1.022839601391910 


c /a 


-36.97653243659754 


-36.97653243473214 


c-i/ao 


-34.01491655155776 


-34.01491655142099 


C-2/ao 


-11.33333333333512 


-11.33333333333343 



5.9 pp ttbb + X 

This process has been compared with the results given 
in [3S] . We have set up the process both in the 't Hooft 
Vcltman scheme and in dimensional reduction and suc- 
cessfully compared the fully renormalised results as an 
internal consistency check. The results below are given 
in the 't Hooft Vcltman scheme, and only the counter- 
terms for \M\1 t Sm are included. 



parameters 





500.0 


N f 


5 


(I 




N f ,h 


1 




174.0 


m b 


0.0 


r t 


0.0 


9s 


1.0 



Using the above parameters and the phase space 
point of Tab. Owe obtain the following results. 



result uu — > ttbb 





GoSam 


Ref. [31] 


a ■ 10 s 


2.201164677187755 


2.201164677187727 


c /a 


8.880263116574282 


8.880263117410131 


c-i/ao 


-4.730495922109534 


-4.730495921691266 


C-2/ao 


-5.333333333333468 


-5.333333333333190 
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E 


Px 


Py 


Pz 


u/g 


250.0 


0.0 


0.0 


250.0 


u/g 


250.0 


0.0 


0.0 


-250.0 


b 


147.5321146846735 


24.97040523056789 


-18.43157602837212 


144.2306511496888 


b 


108.7035966213640 


103.2557390255471 


-0.5484684659584054 


33.97680766420219 


b 


194.0630765341365 


-79.89596300367462 


7.485866671764871 


-176.6948628845280 


b 


49.70121215982584 


-48.33018125244035 


11.49417782256567 


-1.512595929362970 



Table 8 Kinematic point used in pp — ► bbbb. 





E 


Px 


Py 


Pz 


«/s 


250.0 


0.0 


0.0 


250.0 


u/g 


250.0 


0.0 


0.0 


-250.0 


t 


190.1845561691092 


12.99421901255723 


-9.591511769543683 


75.05543670827210 


i 


182.9642163285034 


53.73271578143694 


-0.2854146459513714 


17.68101382654795 


b 


100.9874727883170 


-41.57664370692741 


3.895531135098977 


-91.94931862397770 


b 


25.86375471407044 


-25.15029108706678 


5.981395280396083 


-0.7871319108423604 



Table 9 Kinematic point used in pp — >• ttbb. 



result gg — > ttbb 

GoSam Ref. [38] 

a ■ 10 8 8.279470201927135 8.279470201927128 

co/a 21.83922035777929 21.83922035648926 

c_i/a -12.59181277770347 -12.59181277853837 

c- 2 /a -8.666666666666764 -8.666666666666549 



result gg -> W+W~bb 

GoSam Ref. [39] 

a ■ 10 8 1.549796787502985 1.549795815702494 

c /ao -17.80558461276584 -17.80558440908488 

c_i/a -19.61125131175888 -19.611251301307803 

c_ 2 /a -8.666666666666668 -8.66666666666661 



On an Intel Core i7 950 at 3 GHz the evaluation of a 
single phase space point took 393 ms in the uu channel 
and 12.3 s in the gg channel. The code was compiled 
with gf ortran without optimisations. 



result uu — > W + W bb 

GoSam Ref. [39] 

a ■ 10 s 2.338048681706755 2.338048676370483 

co/a -5.936151367348438 -5.936151368788066 

c_i/a -10.44868110371249 -10.44868110378090 

e-2/ao -5.333333333333312 -5.333333333333336 



5.10 pp^W+W'bb 

The subprocesses uu -> W + W~bb and gg — > W + W~bb 
have been calculated both in [38] and (39]. Accord- 
ingly, the results below are given in the 't Hooft Vcltman 
scheme, where only the counterterms for |A1|^ t Sm are 
included. 



parameters 





500.0 


Nf 


5 


/' 






i 


mi 


174.0 


m b 





r t 





9s 


1 


M z 


91.188 




2.44140351 


M w 


80.419 


r w 





1/a 


132.50686625 







With the above parameters and the kinematics dchncd 
in Tab. [TUlwe obtain the following results. 



5.11 ud — > W + ggg 

The amplitude ud — > W + ggg is an important channel 
in the calculation of the process pp — > W + + 3 jets. The 
QCD corrections to this process have been presented in 
Refs. [1E11]. 

The subprocess with one quark pair and three glu- 
ons consists of more than 1500 Feynman diagrams. We 
have computed the amplitude including the leptonic de- 
cay of the VF-boson and compared our result to |38j . 



parameters 





500.0 


Nf 


5 


/' 




N f,h 


1 


mt 


174.0 


M z 


91.188 


rt, r w , r z 


0.0 


M w 


80.419 


9s 


1.0 


Gf 


1.16639 • 10~ 5 



Furthermore, the values for the dependent parameters 
are cos 2 9 W = /Mf and a = ^/2G F sin 2 6 W /vr. 
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E 


Px 


Pv 


Pz 


u/g 


250.0 


0.0 


0.0 


250.0 


u/g 


250.0 


0.0 


0.0 


-250.0 


w+ 


154.8819879118765 


22.40377113462118 


-16.53704884550758 


129.4056091248114 


w~ 


126.4095336206695 


92.64238702192333 


-0.4920930146078141 


30.48443210132545 


b 


174.1159068988160 


-71.68369328357026 


6.716416578342183 


-158.5329205583824 


b 


44.59257156863792 


-43.36246487297426 


10.31272528177322 


-1.357120667754454 



Table 10 Kinematic point used in pp — ► W+W bb. 





E 


Px 


Py 




u 


250.0 


0.0 


0.0 


250.0 


d 


250.0 


0.0 


0.0 


-250.0 


W+ 


162.5391101447744 


23.90724239064912 


-17.64681636854432 


138.0897548661186 


g 


104.0753327455388 


98.85942812363483 


-0.5251163702879512 


32.53017998659339 


g 


185.8004692730082 


-76.49423931754684 


7.167141557113385 


-169.1717405928078 


g 


47.58508783667868 


-46.27243119673712 


11.00479118171890 


-1.448194259904179 



Table 11 Kinematic point used in ud — > W + ggg. 



For the phase space point of Tab. [TT] we obtain the 
numbers below. 



result ud — > W+ggg 





GoSam 


Ref. [38] 


a ■ 10 7 


8.552735739069321 




c /a 


-36.45372625230239 


-36.4536949986367 


c-i/ao 


-34.70010131004584 


-34.70007155977844 


c-2/ao 


-11.66666666666747 


-11.666656664302845 



On an Intel Core 2 i5 Laptop at 2.0 GHz the evaluation 
of a single phase space point took about 2.5 s for ud — > 
e + v e ggg and about 7.5 s for on-shell W's without decay. 
The code was compiled with gf ortran -02. 

5.12 ud — > W + (— > v e e + )bb (massive 6-quark) 

The process ud — > W + bb, with an on-shell VK-boson, 
has been studied in [100) . while the effects of the W- 
decay have been recently accounted for in jlOlj , and im- 
plemented within MCFM. We consider the latter pro- 
cess, and compare the renormalised amplitude evalu- 
ated by MCFM. The 6-quark is treated as massive in 
all diagrams except in the vacuum-polarisation like con- 
tributions. 



parameters 



/' 


80.0 


9s 


1 


mt 


172.5 


m b 


4.75 


M z 


91.1876 


M w 


80.419 


r w 


2.1054 


Gf 


0.0000116639 


v ud 


0.975 







Using the above parameters and the kinematics given 
in Tab. [T2l we obtain the following results. 





result ud — > v e e 


+ bb 




GoSam 


MCFM-6.0 


a ■ 10 7 


1.884434667673654 


1.88443466774536441 


co/ao 


41.21712989438873 


41.217129894410029 


c-1/00 


26.60367070701196 




C-2/Q.0 


-2.666666666666624 




IR-i 


26.60367070701218 




IR-2 


-2.666666666666667 





The evaluation of a single phase space point took 
9.12 ms on a 2 GHz processor. The code was compiled 
with gf ortran -02. 



6 Conclusions 

We have presented the program package GoSam which 
produces, in a fully automated way, the code required 
to perform the evaluation of one-loop matrix elements 
for mult i- particle processes. The program is publicly 
available at http://projects.hepforge.org/gosam/ 
and can be used to calculate one-loop amplitudes within 
QCD, electroweak theory, or other models which can be 
imported via an interface to LanHEP and UFO, also 
included in the release. Monte Carlo programs for the 
real radiation can be easily linked through the BLHA 
interface. 

GoSam is extremely flexible, allowing for both uni- 
tarity-based reduction at integrand level and traditional 
tensor reduction, or even for a combination of the two 
approaches when required. The amplitudes are gener- 
ated in terms of Feynman diagrams within the dimen- 
sional rcgularisation scheme, and optionally the calcu- 
lation can be carried out cither in the 'tHooft Velt- 
man or in the dimensional reduction variant. The user 
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76.084349979114506 
1998.0331337409114 
-953.55303294091811 
-190.20402007017753 
-417.39085287123652 
-360.80087787946474 



Px 

0.0 
0.0 

955.01676368653477 
194.22279012023398 
468.23544715890415 
456.64248275435313 



Py 

0.0 
0.0 

50.025808060592873 
4.3588877692445251 
208.22173996408185 
-262.60643579391922 



Pz 

76.084349979114506 

-1998.0331337409114 
17.060211586132972 
39.063065018596490 
40.625785184424117 

-96.749061789153586 



Table 12 Kinematic point used in ud — > W+bb. 



can choose among different libraries for the master in- 
tegrals, and the setup is such that other libraries can 
be linked easily. 

The calculation of the rational terms is very modu- 
lar and can proceed either along with the same numer- 
ical reduction as the rest of the amplitude, or indepen- 
dently, before any reduction, by using analytic informa- 
tion on the integrals which can potentially give rise to 
a rational part. In the current version of the code, UV- 
renormalisation counterterms are provided for QCD cor- 
rections only. Further improvements concerning the full 
automatisation of electroweak corrections are planned. 

Different systems to detect and rescue numerical in- 
stabilities are implemented, and the user can switch 
between them without having to rc-gcncratc the source 
code. Due to a careful organisation of the calculation 
both at the code generation stage and at the reduc- 
tion stage, the runtimes for multi-particle amplitudes 
are very satisfactory. Moreover, the GoSam generator 
can also produce codes for processes that include inter- 
mediate states with complex masses. 

Within the context of the automated matching of 
Monte Carlo programs to NLO virtual amplitudes, Go- 
Sam can be used as a module to produce differential 
cross sections for multi-particle processes which can be 
compared directly to experiment. Therefore we believe 
that GoSam can contribute to the goal of using NLO 
tools as a standard framework for the LHC data anal- 
ysis at the TeV scale. 



by the German Federal Ministry of Education and Research. 
The work of G.O. was supported in part by the National 
Science Foundation under Grant PHY-0855489 and PHY- 
1068550. The research of F.T. is supported by Marie-Curie- 
IEF, project: "SAMURAI- Apps" . We also acknowledge the 
support of the Research Executive Agency (REA) of the Eu- 
ropean Union under the Grant Agreement number PITN-GA- 
2010-264564 (LHCPhenoNet). 

Appendix 

A: Examples included in the release 

In the following we give results for the processes listed in 
the examples directory. Unless stated otherwise, we as- 
sume that the coupling constants (e and g s in the stan- 
dard model) have been set to one in the input card. The 
conventions for the returned numbers (ao, cq, c_i, c_2) 
are as stated in Section 12.51 Dimensionful parameters 
are understood to be in powers of GeV. 

As an illustration of the potential of GoSam, we 
display in Table[2]thc timings required by a wide list of 
benchmark processes. The first value provided in the ta- 
ble is the time required for the code generation ( Gener- 
ation, given in seconds) : we remind the reader that this 
operation only needs to be performed once per process. 
The second value is the timing for the full calculation 
of the amplitude at one phase-space point {Evaluation, 
in milliseconds). Results are obtained with an Intel(R) 
Core(TM) i7 CPU 950 @ 3.07GHz. 
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A.l: How to run the examples 

The example directories only define the system inde- 
pendent part of the setup. All settings which are de- 
fined in the file system. rc (see Section [5]) must be 
put either in a file called $H0ME/ . gosam or in the file 
setup, in in the GoSam examples/ directory. A script 
runtests . sh is provided to generate, compile and run 
the test programs. The names of the directories respec- 
tively examples to be run should be specified at the 
command line, e.g. 

. /runtests . sh eeuu bghb 
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parameters 

E 74.7646520969852 /x 2 AE 2 

cj) 2.46 8 1.35 

M z 91.1876 r z 2.4952 

Mvy cos6 w M z sm0 w 0.47303762 



result e+e — > uu 

GoSam analytic 

a 3.7878306213027528 

c /a 1.86960440108932 x C F (tt 2 - 8) x C F 

c_i/o -3.0000000000000 x C* F -3 x C F 

c- 2 /a -2.0000000000000 x C F -2 x C F 



Process 


Generation [s] 


Evaluation [ms] 


bg -> 


236 


2.49 


dd — > tt 


341 


4.71 


dd^ ft (DRED) 


324 


4.05 


dg -» dg 


398 


3.08 


d 3 -> dg (DRED) 


402 


3.28 


e+e - — > tt 


221 


1.27 


e+e" -» if (LanHEP) 


180 


1.27 


e+e - — > mi 


122 


0.65 


33 -> 39 


525 


1.69 


99 -> 93 (DRED) 


428 


1.66 


93 -* 33 (LanHep) 


1022 


1.70 


33 -> 3^ 


529 


15.18 


99 -> tt 


1132 


24.65 


gg -> (DRED) 


957 


30.13 


33 — > tt (UFO) 


1225 


29.45 


H — > 77 


140 


0.24 


9& — > e~i> e t 


337 


2.89 


ud — ¥ e~v e 


71 


0.09 


ud —> e~u e g 


154 


1.15 


uu — > (id 


186 


2.06 


■ud -> W+lV+cs 


1295 


17.37 


77 — > 77 


597 


6.08 



Table 13 Time required for code generation and calculation 
of one phase-space point. The results were obtained with an 
Intel (R) Core(TM) i7 CPU 950 @ 3.07GHz. The time for the 
evaluation of a phase space point is taken as the average of 
the time obtained from the evaluation of 100 random points 
generated using RAMBO [102j. where the code was compiled 
using gf ortran without any optimisation options. The gener- 
ation of the i?2 term was set to explicit. 



If the script is invoked without arguments it will loop 
over all subdirectories. A second script, summarize . sh, 
can be used in order to collect the test results and print 
a summary to the screen. The command 

. / summarize . sh 

will produce an output like the following one. 

+ bghb (succeeded) 
+ eeuu (succeeded) 

grep : ./ddtt/. . . : No such file . . . 

The examples e + e~ — > tt have an explicit dependence 
on the Golem95C library and will therefore fail if the 
extension golem95 is not added. 

A. 2: e + e~ — > uu 

The following parameters and momenta have been used 
to produce the numerical result: 





E 




Px 




Py 




Pz 


e+ 


E 














E 


e~ 


E 














-E 


u 


E 


E 


sin 9 sin </> 


E 


sin 9 


cos <j> 


E cos9 


u 


E 


-E 


sin 9 sin <f> 


-E 


sin 9 


COS <j> 


-E cos e 



A.3: e+e" tt 

This example has been produced twice: once with the 
default model file and once with a model file imported 
from LanHEP [55] ■ Thus it also can serve as an exam- 
ple of how to import model parameters from LanHEP. 
The result is given in dimensional reduction, and no 
rcnormalisation terms are included. 



parameters 



M z 


91.1876 


r z 


2.4952 


M w 


cos W Mz 


sin 6 W 


0.47303762 


m t 


172.5 


V 2 





The following results are obtained with the above pa- 
rameters and the kinematic point of Tab.[T4l 





result e+e - 


-> tt 




GoSam 


analytic 


do 


6.3620691850584166 


6.3620691850631061 


c Q /a 


13.182472828297422 


13.182472828302023 


c-i/ao 


12.211527682024421 


12.211527682032367 


C-2/ao 


0. 


0. 



A. 4: uu -> dd 

This example has been produced twice: once in the 
'tHooft Veltman (HV) scheme and once with dimen- 
sional reduction (DRED). Only the result in the HV 
scheme will be listed below, for the DRED calculation 
see the directory uudd_dred. 

parameters 
fi 91.188 I N f 2 

Using the above parameters and the phase space point 
of Tab.fTSl we obtain the following numbers. 
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E 


Px 


Py 


Pz 


e+ 


74.7646520969852 


0. 


0. 


74.7646520969852 


e~ 


6067.88254935176 


0. 


0. 


-6067.88254935176 


t 


5867.13826404309 


16.7946967430656 


169.437140279981 


-5862.12966020487 


t 


275.508937405653 


-16.7946967430656 


-169.437140279981 


-130.988237049907 



Table 14 Kinematic point used in e+e — > tt 





E 


Px 


Py 


Pz 


u 


102.6289752320661 








102.6289752320661 


u 


102.6289752320661 








-102.6289752320661 


d 


102.6289752320661 


-85.98802977488269 


-12.11018104528534 


54.70017191625945 


d 


102.6289752320661 


85.98802977488269 


12.11018104528534 


-54.70017191625945 



Table 15 Kinematic point used in uu — > dd. 



result uu — > dd 





GoSam(HV) 


Ref. |T03l 


ao 


0.28535063700913421 


0.28535063700913416 


c /a 


-2.7940629929270155 


-2.7940629929268876 


c-x/aa 


-6.4881359148866604 


-6.4881359148866391 


c-2/ao 


-5.3333333333333 


-5.3333333333333 



A. 5: gg -> gg 

This example has been produced both with the default 
model file and with a model file imported from LanHEP. 
Further, it has been calculated in the 't Hooft Veltman 
scheme and in the dimensional reduction scheme. Only 
the results in the 't Hooft Veltman scheme are listed 
below, for further details please see the subdirectories 
gggg_dred and gggg_lhep. The result is for the hclicity 
configuration g(+)g(+) —> g[—)g(—), and pure Yang- 
Mills theory, i.e. fermion loops are not included. 



parameters 



M 2 


442 


N f 





a s 


0.13 







Evaluating the amplitude for above parameters and the 
phase space point given in Tab.lTSlwc obtain the follow- 
ing results. 





result gg ~~, 


99 




GoSam(HV) 


Ref. [1041 




14.120983050796795 


14.120983050796804 


c /a 


-124.0247557942351 


-124.02475579423495 


c-i/ao 


55.003597347101078 


55.003597347101035 


c-2/ao 


-12.00000000000000 


-12. 



A.6: gg -t gZ 

As this process has no tree level amplitude, the result 
is for the one-loop amplitude squared. 



parameters 



M 2 


S12 




1 


M z 


91.1876 


r z 





sin 8 W 


0.4808222 


M w 


cos 6 W Mz 


N f 


2 







With the above parameters and the kinematics given 
in Tab.fTTl we obtain the following result. 



result gg —> gZ 





GoSam 


Ref. [Tu5l 


ao 








0.1075742599502829 


0.10757425995048300 



A.7: dd -> tt 

This example has been calculated in the 't Hooft Velt- 
man scheme and in the dimensional reduction scheme. 
Only the results in the 't Hooft Veltman scheme are 
listed below, for the renormalised amplitude with Nf = 
5 and the top mass renormalised on-shcll. 

For further details please see the subdirectories ddtt 
and ddtt_dred. 



parameters 



m t 


172.5 


M 2 


mf 


a s 


1 


N f 


5 



With the above parameters and the kinematics given 
in Tab.fTSl we obtain the following results. 



result dd — > tt 




GoSam(HV) 


Ref. I27U106I (MCFM) 


a 


0.43024349783870747 


0.43024349783867882 


c /a 


-22.526901042662193 


-22.526901042658068 


c-i/ao 


10.579577611830414 


10.579577611830567 


c-2/ao 


-2.6666666666666599 


-2.666666666666721 
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E 


Px 


Py 


Pz 


Pi 


220.9501779577791 








220.9501779577791 


P2 


220.9501779577791 








-220.9501779577791 


PS 


220.9501779577791 


119.9098300357375 


183.0492135511419 


-30.55485589367430 


Pi 


220.9501779577791 


-119.9098300357375 


-183.0492135511419 


30.55485589367430 



Table 16 Kinematic point used in gg — > gg. 





E 


Px 


Py 


Pz 


9 


100 








100 


9 


100 








-100 


9 


79.2120540156 


3.65874234516586 


- 25.1245942606679 


75.0327786308013 


Z 


120.7879459844 


-3.65874234516586 


25.1245942606679 


- 75.0327786308013 



Table 17 Kinematic point used in gg — > gZ. 





E 


Px 


Py 


Pz 


d 


74.7646520969852 








74.7646520969852 


d 


6067.88254935176 








-6067.88254935176 


t 


5867.13826404309 


16.7946967430656 


169.437140279981 


-5862.12966020487 


J 


275.508937405653 


-16.7946967430656 


-169.437140279981 


-130.988237049907 



Table 18 Kinematic point used in dd — > tt. 



A.8: gg 



tt 



The result is for the renormalised amplitude in the HV 
scheme. 



parameters 



nil 


171.2 


n 





N f 


5 


/' 


71.2 



With the above parameters and the kinematics given 
in Tab.fTOl we obtain the following results. 



result gg — > tt 




GoSam(HV) 


Ref. 1270 1061 (MCFM) 


a 


4.5576116986983433 


4.5576116986983424 


c /a 


15.352143751168184 


15.352143750919995 


c-i/ao 


-27.235240992743407 


-27.235240936279297 


c-2/ao 


-6.0 


-6.0 



A.9: bg -> Hb 

For this process the mass of the b-quark is set to zero. 
However, in order to have a coupling between the b- 
quark and the Higgs boson, the following Yukawa cou- 
pling is implemented in the model file: 



£ yuk = Y m lj)Lll)R4>, y i 



Hb 



m b (n) 



parameters 








m b (n) 


2.937956 


m H 


120 


V 


246.2185 


M 


91.188 







With the above parameters and the kinematics given 
in Tab. [20] we obtain the following results. 





result bg — > 


Hb 




GoSam(HV) 


Refs. |TP71I551 


a ■ 10 7 


2.09926265849001642 


2.09926265848997195 


co /ao 


-24.131948141318752 


-24.131948141995107 


c-i/ao 


11.957924609547224 


11.957924605423791 


c-2/ao 


-5.6666666666666643 


-5.6666666666666670 



A.10: H 



11 



The decay width Fjj-^-y-y of this loop induced process 
is known analytically at lowest order. For comparison 
we used the equations including the top loop and the 



bosonic contribution given in |108U109| . The decay width 
can be expressed as 



G F a 2 m H - 
FH ^- IW2V - r(TW ' Tt) 

where t% = m H /(Amf) for i = W,t. 



(A.l) 



parameters 



m H 124.5 


m t 


172.5 


m w 80.398 






result H 






GoSam 




Refs. 108 : 109 


r{rw,Tt) 3.366785118586698 


3.36678512043889 



27 





E 


Px 


Py 


Pz 


g 


137.84795086008967 


0. 


0. 


137.84795086008967 


9 


3161.1731634194916 


0. 


0. 


-3161.1731634194916 


t 


3058.6441209877348 


16.445287185144903 


165.91204201912493 


-3049.2945357402382 


t 


240.37699329184659 


-16.445287185144903 


-165.91204201912493 


25.969323180836145 



Table 19 Kinematic point used in gg — > it 





E 


Px 


Py 


Pz 


b 


250 








250 


9 


250 








-250 


H 


264.4 


-83.84841332241601 


-86.85350630148753 


-202.3197272300720 


b 


235.6 


83.84841332241601 


86.85350630148753 


202.3197272300720 



Table 20 Kinematic point used in bg — ► Hb. 



A. 11: ud — > e v e 

This example has been calculated in the 't Hooft Velt- 
man scheme and in the dimensional reduction scheme. 
Only the results in the 't Hooft Veltman scheme are 
listed below, for the renormalised amplitude. In addi- 
tion to a calculation with the default model file, calcu- 
lations using LanHEP [55] and UFO [M] are also con- 
tained in the examples directory. 

parameters 



v/s 200 | n 91.1876 

With the above parameters and the kinematics given 
in Tab. [21] we obtain the following results. 





result ud — >• e 






GoSam(HV) 


Ref. [35] 


a 


1.4138127601912656 


1.4138127601912673 


c /a 


5.4861229357937624 


5.4861229357937660 


c-i/oo 


0.18879169932851950 


0.18879169932852413 


c-2/ao 


-2.666666666666667 


-2.6666666666666665 



A. 12: ud — > e v e g 

We list the renormalised amplitude in the HV scheme. 



parameters 



M w 


80.398 


Fw 


2.1054 


sin 9 W 


0.4808222 


M z 


Mw 1 cos W 


N f 


5 


V ud 


0.97419 


l< 2 


312 







With the above parameters and the kinematics given 
in Tab. I2"2l we obtain the following results. 



result ud -— > e u e g 





GoSam(HV) 


Ref. 39 


a ■ 10 7 


2.8398509625435832 


2.8398509625435922 


co /ao 


-8.6052919370147745 


-8.6052919368774248 


c-r/ao 


-18.722010655600936 


-18.722010655557121 


C-2/ao 


-5.6666666666666 


-5.66666666666667 



A. 13: g b — > e v e t 

We list the renormalised result in the dimensional re- 
duction scheme. 



parameters 



M w 80.4190 




2.04760 


M z 91.1876 


r z 


2.49520 


m t 171.2 


n 





H 71.2 


e 


0.30794906326863203 



With the above parameters and the kinematics given 
in Tab. I2"31 we obtain the following results. 





result g b — > e 


u e t 




GoSam 


Ref. I27U106I fMCFM) 


a ■ 10 2 


8.52301540675800134 


8.52301540708130106 


c /a 


-79.879718568538991 


-79.879718569273024 


e_i/ao 


26.570185488790770 


26.570185487963364 


c-2/ao 


-4.3333333333333401 


-4.3333333331689596 



A.14: ud-> W+W+cs 

Results are given for the unrenormalised amplitude in 
the dimensional reduction scheme. 

parameters 
fj, 80 I N f 5 

With the above parameters and the kinematics given 
in Tab. [24] we obtain the following results. 





result W+W+ 


c s 




GoSam 


Ref. [201 v3] 


a 






co /ao 


23.3596455167118 


23.35965 


c-i/ao 


13.6255429251954 


13.62554 


c-i/ao 


-5.333333333333 


-5.33333 
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E 


Px 


Py 


Pz 


u 


100 








100 


2 


100 








-100 


e~ 


100 


75.541566535633046 


30.240603423558878 


-58.128974100026611 


V e 


100 


-75.541566535633046 


-30.240603423558878 


58.128974100026611 



Table 21 Kinematic point used in ud — > e u e . 



500 
500 

483.244841094218 
279.253370247231 
237.501788658551 



Px 




-86.3112218694181 
6.62401666401929 
79.6872052053988 



Py 




147.629518147233 
-5.58083951102529 
-142.048678636208 



Pz 

500 
-500 

-451.975082051212 
279.119009435087 
172.856072616124 



Table 22 Kinematic point used in ud ■ 





E 


Px 


Py 


Pz 


a 


1187.7086110647201 








1187.7086110647201 


b 


2897.148136260289 






-2897.148136260289 


e~ 


2293.0435558834492 


629.81047833131981 


258.58120146220904 


-2189.6399870328105 


V e 


509.48956356743611 


144.72113807954338 


19.883362437475 


-488.098411670514 


t 


1282.3236278741238 


-774.53161641086319 


-278.46456389968404 


968.29887350775562 



Table 23 Kinematic point used in gb — > e u e t. 



B: Explicit reduction of R 2 rational terms 

In this Appendix we list all integrals which give rise to 
R 2 terms as we use these expressions in their explicit 
construction. We use the definition 

H 2e d n qq^ ■ ■■q^(fi 2 ) a 



T n/2 



Di-D 



N 



with 



and 



Di = {q + nf -mf 



(ji — Tj) — m 2 — m 



The integrals up to 0(e) are 
1 

Y 

1 



e-I?'°(S) = -i-Sn 



jn.l 
2 

T n,0 



1 



e-ir(S) = l 



e-l 



(S) 



—g^ {Sn + S 12 + S 22 ) 



(B.2) 



(B.3) 

(B.4) 

(B.5) 

(B.6) 
(B.7) 
(B.8) 



rn,l;/x 



^(5) = --(rf+rf +rf) 



1=1 



7-n,l;pip 2 
J 4 



it 2 (S) 



(S) = T2 9 
1 
6 



e-l 



(S) = ^ [9" 9"} 



(B.10) 
(B.ll) 
(B.12) 

(B.13) 

(B.14) 
(B.15) 
(B.16) 



(B.9) 
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